DX3D.L  .  ■  ■  •^C'^TOX  T i^r.nABY 
NAVAL  POSTGli-ADUATS  SCHOOL 
MONTEREY.  CALIFORNIA  93943-8002 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

DIGITALLY  CONTROLLED  "PROGAMMABLE " 

ACTIVE  FILTERS 

by 

Panagiotis  Andresakis 

December  1985 

Thesis 

Advisor:               Sherif  Michael 

Approved  for  public  release;  distribution  is  unlimited 


T224636 


CURITY  CLASSIFICATION  OF  THIS   PAGE 


REPORT  DOCUMENTATION  PAGE 


REPORT  SECURITY  CLASSIFICATION 


lb.  RESTRICTIVE   MARKINGS 


I.  SECURITY  CLASSIFICATION  AUTHORITY 


.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3     DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


I.  NAME  OF  PERFORMING  ORGANIZATION 

(aval  Postgraduate  School 


6b   OFFICE  SYMBOL 
(If  applicable) 

Code    6  2 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


ADDRESS  (City.  State,  and  ZIP  Code) 


[onterey,    CA      93943-5100 


7b.   ADDRESS  (C/ty,  State,  and  ZIP  Code) 

Monterey,  CA   93943-5100 


NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b.  OFFICE  SYMBOL 
(If  aoolicable) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


ADDRESS  (City,  State,  and  ZIP  Code) 


10    SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 
NO 


TASK 
NO 


WORK   UNIT 
ACCESSION   NO 


TITLE  (Include  Security  Classification) 
IIGITALLY    CONTROLLED    "PROGRAMMABLE"    ACTIVE    FILTERS 


PERSONAL  AUTHOR(S) 

lNDRESAKIS,    Panagiotis 


a    TYPE  OF  REPORT 

laster's    Thesis 


13b    TIME  COVERED 
FROM  TO 


14    DATE  OF  REPORT   (Year.  Month,  Day) 

1985  December 


15    PAGE  COUNT 
201 


SUPPLEMENTARY  NOTATION 


COSATI  COOES 


FIELD 


GROUP 


SUB-GROUP 


18   SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

digitally  controlled  analog  filter 


± 


ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  numtier) 

In  this  research  a  general  purpose  digitally  controlled  analog  filter 
.s  presented.   The  novel  design  is  a  cascade  of  second-order  sections  that 
ire  individually  programmed  to  achieve  any  filtering  toplogies.   Two- 
)inary  words  are  used  to  control  the  pole  frequency  i^   p  and  selectivity 
)p  of  each  section  independently.   Each  second-order  section  is  a 
Jeneralized-Immittance  Converter  (QIC)  biquads  which  are  known  for  their 
ligh  stability  and  low  active  and  passive  sensitivity.   CMOS  switches  are 
ised  to  electronically  relocate  the  minimum  number  of  passive  elements  to 
ichieve  function  programmability .  Switches  are  also  used  to  select  the 
lumber  of  cascaded  sections  to  realize  higher  order  transfer  functions. 


I.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
0  UNCLASSIFIED/UNLIMITED       D  SAME  AS  RPT.  D  OTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


NAME  OF  RESPONSIBLE  INDIVIDUAL 

Sherif  Michael 


22b  TELEPHONE  (Include  Area  Code) 

(408)     646-2252 


22c.  OFFICE   SYMBOL 

Code    62Mi 


)  FORM  1473. 84  MAR 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete. 

1 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Approved  for  public  release;  distribution  is  unlimited 

Digitally   Controlled  "Programmable"  Active  Filters 

by 


Panagiotis  Andresakis 
Lieutanant  Hellenic  Navy 
B.S.,  Hellenic  Naval  Academy  1977 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1985 


"'    """"^     ZTohn  N/.  Dyer 
Dean  of  Science  and  Engineering 


ABSTRACT 
In  this  research  a  general  purpose  digitally  controlled 
analog  filter  is  presented.   The  novel  design  is  a  cascade  of 
second  order  sections  that  are  individually  programmed  to 
achieve  any  filtering  topologies.   Two-binary  words  are  used 
to  control  the  pole  frequency  co  p  and  selectivity  Qp  of  each 
section  independently.   Each  second-order  section 
is  a  Generalized-Immittance  Converter  (GIC)  biquads  which  are 
known  for  their  high  stability  and  low  active  and  passive 
sensitivity.   CMOS  switches  are  used  to  electronically 
relocate  the  minimum  number  of  passive  elements  to  achieve 
function  programmability.   Switches  are  also  used  to  select 
the  number  of  cascaded  sections  to  realize  higher  order 
transfer  functions. 


c' 


TABLE  OF  CONTENTS 


I.   INTRODUCTION 14 

A.  THE  NEED  FOR  AN  ACTIVE  PROGRAMMABLE  FILTER 14 

B.  LINEAR  SYSTEMS 16 

C.  FILTERS  AS  A  SPECIAL  CLASS  OF  LINEAR  SYSTEMS.. 18 

D.  ACTIVE  FILTERS  FUNDAMENTALS 18 

E.  GENERALIZED  IMMITTANCE  CONVERTER 19 

II .   THEORETICAL  ANALYSIS 21 

A.  BIQUADRATIC  TRANSFER  FUNCTIONS 21 

1 .  Low  Pass 22 

2 .  High  Pass 23 

3  .    Band  Pass 23 

4  .    Notch 23 

5.    All  Pass 24 

B .  SENSITIVITY  FUNCTIONS 24 

1 .  Magnitude  Function  Deviation 25 

2 .  Classical  Sensitivity 26 

3  .   Gain  Sensitivity  Product 27 

4.    Determining  the  Variability  of  the 

Transfer  Function  Magnitude 28 


NAVAL  po.-^p':/'-^'^-^^?^"^ 

C.  PROPOSED  GIC  FILTER,  ANALYSIS 29 

1.  Low  Pass  Realization 34 

2.  High  Pass  Realization 36 

3 .  Band  Pass  Realization 37 

4.  Notch  Realization 38 

5 .  All  Pass  Realization 38 

D.  SENSITIVITY  ANALYSIS 39 

E.  STABILITY 52 

III.   PROGRAMMABLE  GENERALIZED  IMMITTANCE  CONVERTER 

FILTER 53 

A.  GENERAL 53 

B.  THE  PROPOSED  GIC  PROGRAMMABLE  FILTER 54 

C.  THE  REALIZED  GIC  PROGRAMMABLE  FILTER 61 

IV.   COMPUTER  SIMULATION 69 

A.  INTRODUCTION 69 

B.  STIMULATION  RESPONSE(S) 70 

1 .  Low  Pass  Realization 70 

2 .  High  Pass  Realization 80 

3 .  Band  Pass  Realization 84 

4  .  Notch  Realization 88 

5 .  All  Pass  Realization 101 


■) , 


V.   REALIZATION  OF  THE  PROGRAMMABLE  GIG  FILTER. ......  105 

A.  EXPERIMENTAL  RESULTS 105 

1.  Low  Pass 105 

2.  High  Pass 107 

3  .   Band  Pass 122 

4.  Notch 122 

5.  All  Pass 142 

B.  CONCLUSION 142 

VI .   COMBINING  HIGHER  ORDER  SECTIONS 149 

VII.   APPLICATIONS  OF  THE  PROPOSED  GIC  PROGRAMMABLE 

FILTER  IN  FREQUENCY  HOPPING  SYSTEMS 167 

A.  BACKGROUND 167 

1.  General  Description  of  FH  Systems 167 

2.  Signal  Generation 168 

a .  Direct  Synthesis 168 

b.  Indirect  Synthesis 170 

B.  PROPOSED  USES 173 

VIII.   CONCLUSION 177 

APPENDIX  A:   PROGRAM  USED  TO  STIMULATE  THE  IDEAL 

GIC  PROGRAMMABLE  FILTER 179 

APPEDNIX  B:   PROGRAM  USED  TO  SIMULATE  THE  NONIDEAL 

GIC  PROGRAMMABLE  FILTER 185 

APPENDIX  C:   PAPER  PRESENTED  AT  THE  19TH  ANNUAL 
ASILOMAR  CONFERENCE  ON  CIRCUITS, 
SYSTEMS  AND  COMPONENTS 191 

LIST  OF  REFERENCES 196 

INITIAL  DISTRIBUTION  LIST 200 

6 


LIST  OF  TABLES 

2.1       Elements  Identification  for  Realizing  the  Most 

Commonly  Used  Transfer  Functions 41 

3.1  The  Elements  Identification  for  Different 

Realization  of  the  GIC  Filter 56 

3.2  The  Truth  Table  of  the  Switches  Logic  Used 

to  Select  the  Filtering  Function 59 

3.3  The  Four  Bitwords  that  control  Fp  and  the 

Corresponding  Capacitor 65 

3.4  The  Four  Bitwords  that  control  Rq  and  Qp 66 

4.1  The  Resulting  Qp  Values  for  the  Different 

Introduced  Rq  Values  in  the  Computer 

Simulation 71 

4.2  The  Resulting   p  Values  for  the  Different 

Introduced  G  Values 72 

4.3  The  25  Possible  Combinations  of  Transfer 

Function  Realizations 7  3 

4.4  Data  from  Figs.  (4.3)  and  (4.4)  Indicating  the 

Affect  of  Frequency  Dependency  of  Ai  (i-1,  2)   79 

4.5  Data  Illustrating  the  Ideal  vs.  the  Nonideal 

Responses  of  the  HPF 86 

4.6  Data  Indicating  the  Ideal  vs.  the  Nonideal 

Responses  of  the  BPF  Realization 97 

7 


LIST  OF  FIGURES 

2.1  Active  Filter  Configuration  Using  the  GIC 30 

2.2  The  Circuit  of  2.1  Prepared  for  Nodal  Analysis.  31 

2.3(a)  The  CGIC  Implementation  Using  Op.  Amps 40 

2.3(b)  Symbolic  Representation  of  the  CGIC  with  40 

Created  3G  and  4G 40 

2.3(c)    The  Basic  Conf irguration 40 

3.1  The  Generalized  Immittance  converter  GIC 

Implementation  Using  Op.  Amplifiers 55 

3.2  Schematic  Diagram  of  the  Programmable  GIC 

Filter  Showing  the  Controlled  Noder 57 

3.3  Different  Elements  Realizations  and  the 

Corresponding  Switches  Used  for  Digitally 

Selecting  the  Filtering  Type 58 

3.3(a)     The  CMOS  Logic  Diagram  Used  to  Control  the 
Analog  CMOS  Switches  of  Fig.  (3.4)  and 
Realize  the  Truth  Table  of  Table  (3.2) 60 

3 . 4  The  complete  Circuit  Diagram  of  the  Programmable 

GIC  Filter 60 

3.5  The  Two  Capacitor  Banks  Realizations  for  the 

Programming  of   p 62 

3.6  The  Resistor  Bank  used  to  Realize  Rq  Needed 

for  the  Programming  of  Qp 63 

8 


3.7  The  Constructed  Blacks  of  Capacitor's  Controlling 

the  f p 67 

3.8  The  Constructed  Block  of  Resistance  Controlling 

QP 68 

4.4       Ideal  LPF  Amplitude  Response 75 

4.2  Ideal  LPF  Amplitude  Response 76 

4.3  Ideal  vs.  Nondieal  CPE  Amplitude  Response  for 

Q=2  and  Vareity  of  Frequencies 77 

4.4  Ideal  vs.  Nonideal  LPF  Amplitude  Response  for 

Frequencies  (7.96k  and  30.1k)  and  a  Variety 

of  Q's 78 

4.5  Ideal  HPF  Amplitude  Response  for  a  Variety  of 

Q's 81 

4 . 6  Ideal  HPF  Amplitude  Response 82 

4.7  Nonideal  HPF  Amplitude  Response 83 

4.8  Ideal  vs.  Nonideal  HPF  Amplitude  Response 85 

4.9  Ideal  vs.  Nonideal  HPF  Amplitude  Response 87 

4 .  10      Ideal  BPF  Amplitude  Response 89 

4 .  11      Ideal  BPF  Amplitude  Response 90 

4.12  .  .  Ideal  BPF  Amplitude  Response 91 

4.13  .  Ideal  BPF  Amplitude  Response 92 

4.14  Nonideal  BPF  Amplitude  Response 93 

4.15  Nonideal  BPF  Amplitude  Response 94 

4.16  Ideal  vs.  Nonideal  BPF  Amplitude  Response 95 

4.17  Ideal  vs.  Nonideal  BPF  Amplitude  Response 96 

4.18  Ideal  Notch  Amplitude  Response 98 

4 .  19      Ideal  Notch  Amplitude  Response 99 

9 


4.20  Ideal  vs.  Nonideal  Amplitude  Response 100 

4.21  Ideal  vs.  Nonideal  Amplitude  Response 103 

4.22  Ideal  vs.  Nonideal  Phase  Response..... 104 

5 . 1  LPF  Amplitude  Response 106 

5.2  LPF  Amplitude  Response.  , 108 

5 .  3  LPF  Amplitude  Response 109 

5 . 4  LPF  Amplitude  Response 110 

5 .  5  LPF  Amplitude  Response Ill 

5 . 6  LPF  Amplitude  Response 112 

5 . 7  LPF  Amplitude  Response 113 

5.8  LPF  Amplitude  and  Phase  Response 114 

5.9  LPF  Amplitude  and  Phase  Response 115 

5.10  HPF  Amplitude  Response 116 

5.11  HPF  Amplitude  Response 117 

5.12  HPF  Amplitude  Response 118 

5.13  HPF  Amplitude  Response 119 

5.14  HPF  Amplitude  Response 120 

5.15  HPF  Amplitude  Response 121 

5  .  16  HPF  Amplitude  and  Phase  Response 123 

5.17  HPF  Amplitude  and  Phase  Response 124 

5.18  HPF  Amplitude  and  Phase  Response 125 

5  .  19  HPF  Amplitude  and  Phase  Response 126 

5 .  20  BPF  Amplitude  Response 127 

5.21  BPF  Amplitude  Response 128 

5 .  22  BPF  Amplitude  Response 129 

5.23  BPF  Amplitude  Response 130 

5.24  BPF  Amplitude  Response 131 

10 


5.25      BPF  Amplitude  Response 132 

5 . 2.6      BPF  Amplitude  Response 133 

5.27  BPF  Amplitude  Response 134 

5.28  BPF  Amplitude  and  Phase  Response 135 

5.29  BPF  Amplitude  and  Phase  Response 136 

5.30  BPF  Amplitude  and  Phase  Response 137 

5.31  BPF  Amplitude  and  Phase  Response 138 

5.32  BPF  Amplitude  and  Phase  Response 139 

5.33  Notch  Filter  Amplitude  Response 140 

5.34  Notch  Filter  Amplitude  Response 141 

5.35  Notch  Filter  Phase  and  Amplitude  Response 143 

5.36  Notch  Filter  Amplitude  Response 144 

5.37  Notch  Filter  Amplitude  Response 145 

5.38  All  Pass  Filter  Amplitude  Response 146 

5.39  All  Pass  Filter  Amplitude  Response 147 

5.40  All  Pass  Filter  Amplitude  and  Phase  Response...  148 

6 . 1  Cascade  GIC  Programmable  Filter 149 

6.2  LPF  Fourth  Order  vs.  Second  Order  Ideal 

Response  from  Computer  Simulation 151 

6.3  Network  with  Logic  Realizing  up  to  the  8th 

Order  LP,  HP,  BP,  N  and  AP  Transfer  Functions  152 

6.4  Fourth  Order  vs.  Second  Order  Experimental 

LPF  Responses 153 

6.5  Ideal  Fourth  vs.  Second  Order  HPF  Responses 

From  Computer  Simulation 154 

6.6  Fourth  Order  vs.  Second  Order  Experimental 

HPF  Responses 155 

11 


6.7  Fourth  Order  vs.  Second  Order  BPF  Amplitude 

Responses  from  Computer  Simulation 157 

6.8  Fourth  Order  vs.  Second  Order  Experimental 

BPF  Responses 158 

6.9  Fourth  Order  vs.  Second  Order  Notch  Filter 

Amplitude  Responses  from  Computer  Simulation.  160 

6.10  All  Pass  Fourth  Order  vs.  Second  Order 

Amplitude  Responses  from  Computer  Simulation.  162 

6 .  11      Chevychev  Response 163 

6.12  Chevychev  Response 164 

6.13  BP-BD  Amplitude  Response 165 

6.14  BP-BD  Amplitude  Response 166 

7.1  Block  Diagram  of  Frequency  Hopping  Modem 169 

7.2  Diagram  of  a  Direct  Frequency  Synthesizer 171 

7.3  Indirect  Signal  Synthesis  Using  a  Phase-Locked 

Loop 172 

7.4  Proposed  Receiver  of  a  F.H.  System  Using  the 

Programmable  GIC  Filter 174 

7.5  Practical  BPF  Realizations  of  the  Composite  GIC  176 


12 


ACKNOWLEDGEMENTS 

Totally  devoted  to  my  father. 

I  also  would  like  to  express  my  gratitute  to  Professor 
Sherif  Michael  for  all  his  teaching  and  understanding. 


13 


I.   INTRODUCTION 

A.   THE  NEED  FOR  AN  ACTIVE  PROGRAMMABLE  FILTER 

The  availability  of  an  analog  filter  with  digitally 
controlled  "programmable"  coefficients  has  been  the  goal  of 
many  researchers  due  to  its  several  attractions.  One 
possibility  of  a  compact,  versatile  analog  filter  under 
remote  control  opens  up  many  novel  and  independent 
application  areas.  Also,  when  a  programmable  filter  is 
combined  with  a  permanent  referenced  memory  which  is 
user-programmable,  this  would  form  an  economical  and 
versatile  device  for  dedicated  stand-alone  applications.  The 
need  for  such  a  device  was  motivated  by  advancement  in  thick 
and  thin  film  technologies  and  continuous  upgrading  of 
systems  specifications  to  take  advantage  of  the  available 
technologies  to  the  limits. 

Linear  analog  filtering  finds  many  applications,  such  as 
speech  processing  (recognition  or  synthesis),  geology, 
instrumentation,  communications,  process  control,  adaptive 
balancing,  etc.  There  has  been  much  emphasis  on  performing 
the  filtering  function  digitally,  largely  because  of  the  ease 
of  varying  and  optimizing  the  transfer  function.  However, 
for  many  reasons,  such  as  cost,  size,  signal  processing 
complexity,  and  bandwidth,  it  would  be  desirable  to  perform 
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the  filter  function  with  linear  components  yet  retain  the 
flexibility  of  varying  the  filter  parameters  digitally. 

Recently,  the  advantages  of  combining  linear  components 
(operational  amplifiers,  resistors  and  capacitors)  and  non 
linear  elements  (switches)  have  been  demonstrated  using 
switched  capacitor  techniques  [36-38].  In  this  research,  we 
are  presenting  the  results  of  realizing  a  continuous  active 
device  using  linear  elements  and  switches  controlled  by 
digital  signals  to  achieve  a  fully  programmable  filter  [32]. 
Several  programming  features  of  the  proposed  filter  are 
reported.  The  first  feature  is  the  ability  of  the  network  to 
realize  the  most  common  filtering  functions  (function 
programmability)  namely:  Low  Pass  (LP),  Band  Pass  (BP),  High 
Pass  (HP),  All  Pass  (AP)  and  Notch  (N)  functions,  using  the 
minimal  set  of  elements.  The  second  feature  is  the  ability 
of  the  network  to  program  (independently)  the  key  parameters 
of  the  filtering  function  chosen  (parameter  programmability) 
namely:  the  pole  resonant  frequency  {o^p)  and  selectivity 
(Qp).  Finally,  the  ability  to  program  the  network  to  cascade 
several  sections  to  achieve  higher  order  filer.  All  of  the 
above  programmability  features  are  performed  independently  to 
realize  a  universal  filtering  network.  In  order  to 
demonstrate  the  idea  of  this  research,  it  is  necessary  to 
introduce  some  theoretical  back  ground. 
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B.   LINEAR  SYSTEMS 

The  box  at  Fig.   (1.1)  illustrates  the  concept  of  the 

linear   system.    In   the   time   domain^   the   system  is 

characterized  by  its  impulse  response  y(t),  which  is  the 

output  signal  y(t)  produced  in  response  to  an  impulse.  for 

an  arbitrary  input  signal  u(t),  the  output  signal  y(t)  is 
given  by  the  well  known  convolution  integral. 


^LL^r  caLL-1)    IlLCi 


C^  \) 


■U.LL) 

^Ll) 

Uls^ 

ic.. 

Fig.  1.1  -  Black  Box  Concept  of  Linear  Syst 


em 


The   system   transfer   function   G(s)   is   the   Laplace 
transform  of  g(t),  thus: 


(sL'.\  ^  ^  Li)     (ixp  C-t.  L^ 


(1.2) 


iwhere.    S=ir  +  iu) 
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Under  the  laplace  transformation  1.1)  becomes: 


ic,^~.G.L'.^-\]L5^  (1-3) 


so  that  the  transfer  function  G(s)  is  the  ratio  of  the  output 
variable  to  the  input  variable. 


^L.^.ji^^  (1.4) 


The  most  general  types  of  linear  systems  are  consisting 
of  a  finite  number  of  lumped,  linear  and  time-invariant 
elements.  The  system  is  characterized  by  an  Nth  order 
ordinary  linear  differential  equation  which  results,  in  most 
general  cases,  a  transfer  function  G(s)  which  is  a  real 
rational  polynomial  function  of  the  complex  variable  S.  Thus 
we  can  write  G(s)  in  general  as 


(s 


c«  = 


^^^)  .  ?hS'\?.-.s"t  -  -.fo   ,^.£^S 


tcs)      t,s^.U1sn..  -vto       ±(<b' 


(1.5) 


K-O 


where  Pk  and  Ik  are  real  numbers  so  that  G(s)  is  real  for 

real  s,  and  the  roots  of  the  polynomial  P(s)  and  E(s)  must  be 
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real  or  occur  in  conjugate  pair. 

Also  by  proper  multiplications  G(s)  can  take  the  form 


p        TrCS-53K) 
6^s..lf!L.^  .^^ (1.6) 


where  C  is  a  constant  extracted  from  G(s)  such  that  E(s)  and 
P(s)  become  monic  polynomials  (leading  coefficients  equal  to 
unity)  ,  Szk  are  the  transmission  zeros  and  Snk  are  the 
natural  modes  at  the  system. 

C.  FILTERS  AS  A  SPECIAL  CLASS  OF  LINEAR  SYSTEMS 

Linear  systems  can  be  distinguished  into  "SPECTRAL 
SHAPING  NETWORKS"  and  "FILTERS."  The  role  of  filters  is  one 
of  selecting  signals  while  the  role  of  spectral  shaping 
networks  is  that  of  modifying  the  input  signal  spectrum  in  an 
arbitrary,  but  predescribed  manner.  Specifically,  we  desire 
that  a  filter  should  do  as  little  as  possible  shaping  on 
signals  in  its  passband;  any  shaping  it  is  considered  a 
distortion  of  the  signal.  On  the  other  hand,  networks  which 
perform  pulse  forming  fall  within  the  spectral  shaping 
category. 

D.  ACTIVE  FILTER  FUNDAMENTALS 

The  distinction  between  passive  and  active  filters  is 
that  the  first  do  not  require  a  power  source  to  perform  their 
function  while  the  second  do.   The  motivation  behind  active 


RC  filters  lies  in  the  desire  to  have  inductorless  filter 
realizations.  It  is  well  known  that  of  the  three  passive 
R,C&L  elements,  the  inductor  is  the  most  nonideal  one.  This 
is  especially  true  at  low  frequencies,  where  inductors  become 
quite  bully  and  have  increased  losses  or  equivalently  lower 
Q-f actors. 
E.   GENERALIZED  IMMITTANCE  CONVERTOR 

One  of  the  methods  of  active  RC  filters  design  consists 
of  simulating  the  inductances  in  the  LC  ladder  by  active  RC 
networks.  This  simulation  can  be  based  on  the  principle  that 
we  want  to  find  a  one  port  network  having  an  input  impedance. 

Zll  =  s*L 

Various  active  elements  as  well  as  synthesis  procedures 
employing  them  have  been  proposed  [1-9].  A  partial  list 
includes : 

(1)  Negative  Impedance  Converter  (NIC) 

(2)  Negative  Impedance  Inverter   (NIV) 

(3)  Postive  Impedance  Converter   (PIC) 

(4)  Gyrator 

(5)  Generalized  Impedance  Converter  (GIC) 

( 6 )  Curent  Conveyor  ( CC ) , and 

(7)  Operational  Amplifier  (OA) 

Although   the   introduction   of   these   elements   has 

stimulated  research  in  the  area  of  active  network  theory, 

very  few  elements  have  made  their  way  to   large   scale 
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production  to  become  available  as  off-the-shelf  items.  The 
reason  for  this  is  mostly  an  economic  one.  For  a  device  to 
become  available  at  low  cost,  it  has  to  be  used  in 
substantially  large  quantities.  It  follows  that  such  a 
device  has  to  be  versatile  enough  to  be  of  use  in  a  number  of 
applications,  of  which  active  filter  design  is  only  one.  The 
analog  circuit  design  area  has  found  these  attributes  in  the 
IC  operational  amplifiers  (Op.  Amp. ) 

The  IC  Op.  Amp.  is  currently  the  most  popular  linear 
active  element.  It  is  available  from  a  large  number  of 
manufacturers,  at  reasonable  cost  and  with  good  performance 
characteristics.  Furthermore,  elements  engineers  have  become 
accustomed  to  the  use  of  Op.  Amp.  It  is  therefore,  only 
natural  that  the  Op.  Amp.  is  becoming  the  most  popular  active 
element  in  the  design  of  active  RC  filters,  and  can  be  found 
in  NIC's,  PIC's,  QIC's,  and  other  circuit  realizations. 
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II.   THEORETICAL  ANALYSIS 

A.   BIQUADRATIC  TRANSFER  FUNCTIONS 

The  filter  as  a  special  class  of  linear  system  has  a 
transfer  function  expressed  in  a  polynomial  quotient  form 
given  as 


p^s)     r,v,  s^t";.,  s'l     +?c 


'JA   t     c  P  P 


fK  :> 


Cb^-   ■  -^  -_:: _r^ 1  ._  K^ 


(2.1) 


where  Pk  and  gk  are  real  numbers  to  that  T(s)  is  real  for  s, 
and  the  roots  of  the  polynomials  P(s)  and  E(s)  must  either  be 
real  or  occur  in  conjugate  pairs.  Also,  in  general,  the 
degree  of  the  numerator  (deg.[P(s)]  =  M)  is  less  than  or 
equal  to  the  degree  of  denominator  (deg[E(s)]  =  N)  and  the 
roots  of  E(s)  are  in  the  open-half  S-plane.  The  E(s)  is 
known  as  the  characteristic  polynomial  or  natural  mode 
polynomial  of  the  linear  system,  and  the  degree  of  E(s),  that 
is  N,  is  the  order  or  degree  of  the  system. 

A  general  second-order  transfer  function  or  "biquad" 
function  may  be  written  as 
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r^^^  ^  P3  s^-^f^Si-Po      _  Pes; 


-^-^lSv%        tc5.;  (2.2) 


where  P(s)  is  the  loss-pole,  or  more  appropriately  here,  the 
transmission-zero  polynomial,  and  E(s)  is  the  natural  pole 
polynomial  mode  as  discussed  above.  It  is  a  usual  practice 
to  express  the  denominator  in  terms  of  _p  and  Qp,  where  _p  is 
the  natural-mode  or  resonance  frequency  and  Qp  is  the 
natural-mode  or  quality  factor.   Thus  (2.2)  becomes 


\  c^\,   _  ^-3  S*^-*-  f/vS  r  \'o 


.^      ^P         ,.z  (2.3) 

The  numerator  coefficients  determine  the  location  of  the 
transmission  zeros  and  hence,  the  type  of  filter  function  the 
biquad  provides.   Special  cases  of  interest  are: 

1.   Low  Pass  (LP) 

For  which  P1=P2=0,  thus  two  transmission  zeros  are  at 
infinity 


Tcs^^       ^o 


Q^    ^"^      ^  (2.4) 
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2.   High  Pass  (HP) 

For  which  pO=pl=0,  thus  two  transmission  zeros  are  at 
infinity,  and 


(Ts)- 


(2.5 


.2       UOp 


5    -f S  V  Ul3  a^ 


3.      Band-Pass    (BP) 

For  which  Po=P2=0,thus  one  transmission  zero  is  at 
infinity  while  the  other  is  at  the  origin,  and 


rc,>-_ ^ 


5  +- L  S  -t-  UO  -T 


(2.6) 


Notch  (N) 


For  which  P1=0  and  the  two  transmission  zeros  are  at 
or  high-pass-notch) ,  leading  to 


'-  joc^        i^a<  UL3p  (depending  if  we  have   low-pass-notch 


\C,~),l 


^-2  .^      2. 

5     +    ^O 


.^^    ^P 


(2.7) 


S-.^S.o^p^ 


2S 


5.   All  Pass  (AP) 

For  which  the  pair  of  zeros  are  at  the  mirror  image 
location  of  the  pair  of  poles,  that  is 


^  (2.8) 


s^,^^  S-^P^ 


^f 


B.   SENSITIVITY  FUNCTIONS 

A  concern  about  the  design  of  a  filter  is  how  close  the 
resulting  response  will  be  to  the  ideal  or  desired  function. 
The  reason  for  response  deviation  from  the  ideal  is  the 
finite  tolerances  of  the  RC  elements,  as  well  as  the  nonideal 
performance  of  the  active  elements.  In  the  latter  case,  not 
only  the  gain  changes  or  tolerances  have  to  be  considered  but 
the  effect  of  the  "limited  amplifier  bandwidth"  on  the  filter 
response  must  also  be  evaluated.  Although  effects  of  initial 
component  tolerances  may  be  "trimmed  out"  during  the  initial 
filter  alignments  or  tuning  process,  a  sensitive  design  will 
deviate  from  the  required  specifications  as  time  process,  due 
to  component  variations  with  temperature,  aging,  humidity, 
etc.  Note  also  that  a  sensitive  design  might  be  extremely 
difficult  to  tune  in  the  first  place,  or  the  initial 
adjustment  will  be  quite  uneconomical. 

The  answer  to  tolerance  question  can  be  obtained  through 

sensitivity  studies.   Considerable  emphasis  has  been  placed, 

in  the  active  filter  literature,  on  the  study  of  sensitivity 
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functions  and  relations.   some  of  the  most  useful  and  widely 
accepted  sensitivity  functions  are  the: 
1.   Magnitude  Function  Deviation 

Assure  a  filter  designed  to  meet  a  certain  magnitude 
characteristics  T(s)  or  T(j'>^).  One  is  concerned  with  the 
deviation  in  |T(jO)j,  that  iSj\.jT(  j\,i^)  j  both  in  passband  and  in 
stopband.  Usually  it  is  desirable  to  express  the  expected 
deviation  in  dB.  The  deviation  D(iO)  dB  in  the  magnitude 
function  may  be  evaluated  as  follows.  Let  the  function  jT(j->^i"; 
change  to  [jT(jai)j  ^     jT(j  )j],  then 


D^a.)  z:  QO  locj  I  ""C^^^^^ATcj^  ^  ^^^^ 


CjU3) 


2.9 


or 


Dlld^  -  8  68  ^ 


1^ 


4.  + 


Ah  (juo) 


U^) 


^       .(^'^) 


(2.10) 


and  for  small  variability  (2.10)  can  be  approximated  as, 


Pluo) 


8.68 


.(^^) 


(2.11) 
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Thus,  the  deviation  in  the  magnitude  response  in  nepers 
is  equal  to  the  per  unit  variability  in  the  magnitude  of  the 
transfer  function.  The  problem  now  reduces  to  that  of 
evaluating  the  per  unit  change  in  jT(jvP)j.  This  is  not  an 
easy  problem  since  T(juj)  is  a  function  of  many  elements  with 
different  tolerances  and  tolerance  statistics.  Furthermore, 
the  per  unit  change  is  function  of  frequency. 

2.   Classical  Sensitivity 

Lets  recall  the  definition  of  the  classical 
sensitivity,  S^  where  y  is  a  variable  of  interest,  usually  a 
function  of  many  parameters  of  which  x  is  one,  then 


->..  ^ —  =  (2.12) 


A 


^>^     '^    V>l^n^) 


Note  that  from  the  above  definition,  S^  is  the  limit  to  as 
Dx  -->0.   thus,  for  small  variations. 


-^^  -   ,   ^  (2.13) 


The  usefulness  of  the  classical  sensitivity  function 
is  evident  from  (2.13).   The  per  unit  or  percentage  change  in 
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y,  due  to  a  given  per  unit  or  percentage  change  in  x,  can  be 
easily  obtained  by  multiplication  with  S,    i.e.. 


h_]^S^   ^  ^ 


"^     J  K   V   A 


(2.14) 


3 .   Gain  Sensitivity  Product 

An  important  consideration  in  the  evaluation  of  the 
sensitivity  of  a  filter  parameter  as  considered  in  [2]  with 
respect  to  the  closed  loop  gain  is  the  tolerance  on  the 
closed  loop  due  to  the  open  loop  gain  variability. 

Thus,  the  gain-sensitivity-product,  G.S^  is  defined 
as: 


(s.   b^  /^  K  S^  (2.15 


where  k  is  the  closed  loop  gain.   We  can  extend  (2.15)  to  the 
open  loop  gain  as: 


(s  S  <^   ^o   S.  ^  "  (2.16) 


and  also  we  can  note  that: 

27-  • 


G.S%65^ 


(2.17) 


and  thus  then  ultimate  good  is  the  variability  or  tolerance 
rather  than  the  sensitivity,  the  gain-sensitivity  product  is 
a  better  index  for  comparing  different  designs. 

4.   Determining  the  Variability  of  the  Transfer  Function 

Amplitude 

Assure  that  the  active  filter  has  t      resistors, 

iTL  capacitors  and  n  amplifiers.   Let  the  amplifier  k  have  an 

open  loop  gain  Aok  and,  possibly,  a  closed  loop  gain  Kk.   The 

variability  of   the  magnitude   function  as  given  by  the 

Reference  [10]  is: 


2     \    ^^       )  ^—        '-\ 


-  gC-'I'IK. 


A  ^ 


(2.18) 


Note  that  each  of  the  sensitivity  functions  is  (2.13) 
is  a  function  of  frequency.  In  a  high  order  filter 
realization,  the  different  sensitivity  functions  might  be 
difficult  to  evaluate. 
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C.   PROPOSED  GIC  FILTER  ANALYSIS 

In  order  to  obtain  the  transfer  functions  of  the  proposed 
programmable  filter  [13]  shown  at  Fig.  (2.1),  nodal  analysis 
was  used  as  follows: 

(1)  The  circuit  of  Fig.  (2.1)  was  replaced  by  the  one  of 
Fig.  (2.2)  in  which  the  two  operational 

cimplifiers  were  replaced  by  the  two  equivalent  dependent 
voltage  sources  as  can  be  seen,  every  element  (node, 
admittance,  voltage  source,  etc.)  was  labeled  and  every 
element  was  associated  to  a  current  direction  and  a  voltage 
polarity. 

(2)  The  kirchoff  current  low  was  written  for  every  node 
except: 

(a)  The  reference, 

(b)  any  node  connected  to  the  reference  by  a  voltage 
source. 

Node  "      K.C.L. 

1  i7=-i3  (2.19) 

2  i2=i5+i6  (2.20) 
3 

4  — 

5  14=17+18  (2.21) 
6 

7 
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Fi.;'.  2.1 


Active  Filter  Configuration  Using  the  GIC 
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]  ■- 


Fiq.  2.2  -  The  Circuit  of  2 . 1  prepared  for  Nodal  Analysis 
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( 3 )  Every    admittance    current  was    expressed    in    terms    of 
nodal  voltages: 

il   =  Yl[V3-Vl]  (2.22) 

i2   =  Y2[V3-l)2]  (2.23) 

i3    =  Y3[a/4-Vl]  (2.24) 

i4   =  Y4[V4-V5]  (2.25) 

i5   =  Y5[T/2-V6]  (2.26) 

i6   =  Y6[V2]  (2.27) 

17  =[Y7[^^5-V6]  (2.28) 

18  =Y8[iJ5]  (2.29) 

(4)  The  source  dependencies  were  listed  expressed  in  terms 
of  nodal  voltages. 


(2.30) 


(2.31) 


and  then  (substitute)  where  necessary. 

(5)  A  matrix  equation  having  the  unknown  voltages  V3 ,  V4 , 
and  V5  related  with  the  desired  transfer  functions  Tl,  T2, 
T3 ,  of  the  filter  was  obtained  as: 


^^-1  ^2 


/A3  ^^ 


0 


^z 


iA 


.^^-^1 


Ai  ^^       At 


-^1'%-M 


k33 


^. 


'\3k 


32 


O 


-\'sU, 


y  ?'0,^ 


2.32) 


V>' 


and  by  substituting  \^z  .2.         T^r  — 

the  above  matrix  equation  takes  the  following  form: 


and     Tb--^-:^ 


I       AtJ       As 


A2     Ai 


0 


.\|2  _xj_  ic 

At       A\       /\\ 


^4 


~^i-<l 


-ib-\^is 


-^1_^8-Ha 


\^ 


T- 


T3 


0 


^\5    I 


■<l 


(2.33) 


So  the  above  matrix  equation  gives  the  3  different 
responses  of  the  programmable  filter  as  functions  of  the  Yl , 
.  .  .  Y7  admittances  and  Al,  A2  (Gains)  the  two  Op.  Amps.  In 
ideal  case  [Al=A2=oo]  the  equation  takes  the  form: 


0   x^ 


T-v 


T- 


^5 


0 


"1. 


-\^5 


[''\ 


(2.34) 


which  expressed  the  different  realizations  discussed  by  [31]. 
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Now  we  can  express  the  above  matrix  equation  for  the 
nonideal  case  in  which  Al  and  A2  are  finite  and  frequency 
dependent  ^>(-\m05  ,  Ki:\)^^l5  where  Wl  and  W2  are  the  Gain 
Bandwidth  Products  (GBWP)  of  the  Op.  Amps.  Al,  A2 , 
respectively. 


HvWr^    i   .^::^   S 


v»J2j        ^: 


0 


X5 


V4 


:is~i.       !|^^^.^i 


•^^V^^-^^ 


U(.^) 


\"3C0 


A-^- 


(2.35) 


1.   Low  Pass  (LP)  Realization 

Can  be  obtained  through  T2(s)  by  substituting 
Y6=Y7=0,  and  if  we  consider  the  ideal  case  where  (Al=A2=oo) , 
then 
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.   /•  /    ;  \  (2.36) 


If  we  further  substitute  the  values  of  the  remaining 
admittances  as  proposed  in  [4],  that  is  Y1=G1,  Y2=sC2, 
Y3=sC3+G3,  Y4=G4,  Y5=G5,  Y8=G8  then  (2.36)  takes  the  form 


(2.37) 


2^SV- 


<^^    ^-^5  C<5-4+G-b^ 


GvGa   G5t-6-3G-fe6iSi-^5'  ^a<^3  S^ 


which  is  the  form  of  a  low  pass  transfer  function  as 
indicated  by  (2.4)  where 
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lo  r    6,    GsHsA^G-H^  I  CzsC^Gq 


(2.38) 


ODp^-    ^'6-466 


C.C.G^  (2.39) 


and 


^?       G 


Qp       c 


1  (2.40) 


p  ^3 

2.   High  Pass  (HP)  Realization 

This  one  is  obtained  through  Tl(s)  substituting  the 
following  values  of  the  admittances:  Y1=G1,  Y2=G2,  Y3=sC3, 
Yr=G4,  Y5=0,  Y6=G6 ,  Y7=sC7,  and  Y8=G8 ,  and  if  we  assume  ideal 
case  (Al=A2->   00),  the 


T1CO-   S  C^C^f-  CC':x^Ga^ 


C%(^\G\ ^C3 G^Q^  S ^GjiC^C  1-  S' 


expresses  the  H.P.  filter  of  (2.5)  where 


(2.41) 


G2C1  C^- 


(2.42) 


(2.43) 
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and 


p  _  (5x^6-4 

&5  _         (2.44) 


3.   Band  Pass  (BP)  Realization 

This  is  also  derived  from  Tl(s)  by  substituting  the 
following  values  of  admittances:  Y1=G1,  Y2=G2 ,  Y3=sC3, 
Y4=G4,  Y5=0,  Y6=G6,  Y7=G7,  and  Y8=sC8 ,  and  if  we  assume  ideal 
case  ( Al=A2-->oo) ,  then 


<s^i  G^  (s6  +  s  (^^6  ^ C':?  i- <C5  Cg  G :i s^ 

J 

which  is  the  form  of  a  BP  transfer  function  as  given  by  (2.5) 
where 


^^  ' 7 (2.46) 

^2  <-^Qv 


^    -  (2.47) 


^?      c 


Oncl 


b 


ri- (2.48 
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4.   Notch  (N)  Realization 

T2(s)  expressed  the  Notch  response  if  the  following 
substitutions  have  been  made.  Y1=G1,  Y2=G2 ,  y3=sC3,  Y4=G4 , 
Y6=0,  Y5=G5,  Y7=sC7,  Y8=G8 ,  and  again  assuming  ideal  case 
(Al=A2-->oo) , 

Xors)-    ^^(^3^1-6:^  ^^>65^^4g^^  _      •  (2.49) 

which  is  the  form  of  a  N  transfer  function  as  defined  by 
(2.7)  where 


>  r    r ^2.50) 


ODp^-  _6\_^;4_G5____  (2.51) 


(2.52 


^?   (Sb 


^9     .   Cj. 


5.   All  Pass  (AP)  Realization 

This   is   derived   from  Tl(s)   with   the   following 
admittances   substitutions.    Y1=G1,   Y2=G2,   Y3=sC3 ,   Y4=G4, 
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Y5=G5,  Y6=0,  Y7=sC7  and  Y8=G8 ,  and  if  once  more  we  assume 
ideal  case,  then 


which  is  the  response  of  an  All  Pass  filter  as   (2.8) 
indicates  where: 


(SO? 

Qp 

=  ~ —     C  ^or 

O\od 

^p^= 

nooiT)) HiiDuo  ^^os-e.  ^  (2.54) 


:> 


(2.55 


Table  2.1  shows  all  the  realizations  proposed  by 
[31].  In  our  research  for  designing  a  programmable  filter 
the  No.  1,  3,  7,  9,  and  12  realizations  were  used  since  they 
offer  the  minimum  admittance  elements  .  change  to  shift  from 
one  to  another. 
D.   SENSITIVITY  ANALYSIS 

Consider  the  CGIC  circuit  shown  in  Figure  (2.2(a))  [12]. 
Assuming  ideal  Op.  Amps.,  the  chain  matrix  of  the  CGIC  can  be 
obtained  as 
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(c)   ^ 


Fig.  2.3Ca)The  CGIC  Implimentation  Using 
Op.  Amps. 

(b)  Symbolic  Representation  of  the 
CGIC  with  created  ports  3G  and  4G 

(c)  The  Basic  Conf irguration. 
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Table  2.1  -  Elements  identification  for 
Realizing  the  Most  Commonly  Used  Transfer  Fuctions 
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L^l  . 


i         0 

(2.56) 

0   K.CS) 


where  h(s),  the  admittance  conversion  function,  is  given  by: 


kc<.)-.  iiM  ]i\  U  (2.57) 


Two  new  ports  can  be  created  across  terminals  3G  and  4G 
as  shown  symbolically  in  Fig.  (2.3  (b)). 

A  synthesis  procedure  is  now  described  which  uses  the 
configuration  of  Fig.  (2.3(c)).  The  transfer  functions 
between  the  input  and  output  terminals  2 ,  3 ,  and  4  are 
readily  obtained  as 

^4U,.  2\7-.lhluUu)_U\-^U^   .U)U\1  Ocs)  (2.59) 


'hJ^)L  ^  T?>  r  Kh-kii)   H>'s!Dc>:)  (2.60) 
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where 

\)c^^  r  <5vN6  ^kcs)Chvi'&^  "       (2.61) 


The  conversion  function  h(s)  and  Y5-Y8  can  be  selected  in 
many  different  ways  and  it  is  found  that  any  second-order 
transfer  function  can  be  realized  [12]. 
Letting 

\I-  ^C-L    ^Gx 

where  i=l,  2,    2,    or  4,  we  have  from  (2.57) 


kc5)^  CsC^^Qz^jCsCz'tQ^)   iCsC.  ^r^O  LSC4+C'4)     (2.62) 


Clearly  by  omitting  one  or  more  conductances  and/or 
capacitances  a  number  of  specific  conversion  functions  can  be 
generated. 

Most  frequently,  filters  are  designed  by  using 
Butterworth,  Chebychev,  Bessel  or  elliptic  approximations  in 
which  the  transmission  zeros  are  located  at  the  origin, 
imaginary  axis  or  at  infinity.   Consequently,  the.  transfer 
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function  can  be  expressed  as  a  product  of  a  second-order 
transfer  functions  of  the  form 


cs^   r  _r±-±_l_ni^Liri^ .  (2.63) 


Oz   S  +Oi  S  ■»-  CTo 


"2 


where  al=a2=0,  ao=al=0 ,  ao=a2=0  or  al=0  for  low  pass  (LP), 
high  pass  (HP),  band  pass  (BP)  or  notch  (N)  section, 
respectively. 

The  coefficients  ai ' s  of  T(s)  for  these  sections  are  all 

positive.   These  sections  can  be  realized  by  choosing  h(s)  in 

2 
a   simple  manner   such  as   k^s,   ^o^'   ^3^^  ''"^4^   °^   their 

reciprocals.  The  different  ki's  ( i=l ,  2,  3,  or  4)  are 
positive  constants.  By  comparing  (2.58-2.61)  with  (2.63), 
circuits  1-10  in  Table  2.1  can  be  obtained.  Circuits  3,  4 
and  7  can  be  regarded  as  realizations  of  simple  RLC  networks 
[31]. 

All  pass  transfer  functions  are  often  needed  for  delay 
equalization  and  these  can  be  realized  by  using  second-order 
transfer  functions  of  the  type  give  by  (2.63)  where  a2=b2, 
al=bl,  and  ao-bo.  Second-order  sections  of  this  class  can  be 
obtained  from  circuits  11  and  12  of  Table  2; 1. 

Figures  (2.3(a)),  (2.3(b)),  and  Table  2.1  show  that  with 
the  exception  of  circuit  10,  the  response  is  obtained  from 
the  output  of  an  operational  amplifier.   Owing  to  the  low 
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output  resistance  of  the  amplifier,  any  number  of  sections 
can  be  cascaded  without  isolating  amplifiers. 

An  important  criterion  of  a  realization  is  its 
sensitivity  to  element  variations.  The  pole  Q  factor  and  the 
undamped  frequency  of  oscillation  from  the  transfer  function 
of  (2.63)  are  defined  as 


Ol^-\bob\  Ibi         ^      '^p-\''ooib2  (2.64) 


For  a  Notch  section,  the  Notch  frequency  is  defined  by 


(2.65) 


and  the  multiplier  constant  can  be  taken  to  be 

HurfDiclbc        br         QilbT.  (2.66) 

for  0C)n>oip  or  ain<obp,    respectively. 

Similarly  for  the  LP,  HP,  and  BP  sections 

HLP.-QoilOo^     HuPraalbz      Q«^  l'\&P  ^Ojt)!  (2.67) 
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For  an  All  Pass  Section,  let 


Q 1 


The   sensitivity   of   a   quantity   x   with   respect   to 
variations  in  an  element  e  is  given  by 


^\  ^1-  "^ 


^   tJ: 


e. 


For  ideal  amplifiers,  the  use  of  (  2  .  64  ) -(  2.  68)  and  Table 
(2.1)  leads  to 


O^lSe!^^ 


e  [^  "  (2.69) 


where  x  represents  any  one  of  the  quantities  defined  by 
( 2.64 ) -( 2 .  68 )  and  e  represents  any  capacitance  of 
conductance.   In  addition,  it  can  be  shown  that 


ZS^o 
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and 


ZlS,^|-2lsf  ]  =  4  ,2.70, 


For  amplifiers  with  a  finite  open-loop  gain  A,  according 
to  [31]  the  circuit  of  (Fig.  1(c))  gives 


^^h.  --  ^^^'.:)o.- 


(2.71) 


where  k=2,  3,  4,  and 

rc>D=V"i  ri  +  f2^3t  [urOfufiUi^lAi^fi/A^  ^\\l^^^2^\i'■\,^2^    1^      (2.72) 

Consider    realizations    in    which    h(s)     =    kls,     such    as    the 
circuits   7(BP),    3 (HP),    9,    10 (N),   where 

(2.73) 


For  real  amplifier  gains  such  that 


1-^1^)  gives  ., 
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(2.74) 


From  (2.64),  (2.73),  and  (2.74),  the  Q-f  actor  and  the 
undamped  frequency  of  oscillation  can  be  obtained  as 


I  ^  C3A^^($,^5  ^. 


C3  A2A0 


\       ' 


=  ?c.= 


1  r 


X5  (^4 


^ 


''  \ 


(2.75) 


r2^o  r 


(2.76) 


where 

The  sensitivities  of  Qpa  and   pa  with  respect  to  the 
amplification  Ao  can  be  written  as 


'K- 


3A( 


'M    ^'K±_     2isA      Q^^y    K5_ 
"m  ^  ^'3   ^^      c^    y^3 


~\ 


(2.77) 


The  use  of  (2.77)  and  (2.78)  leads  to 


Q 


S. 


^'>ib[^-':^y--o^('^^^\^<^^ky'^\ 


(2.78) 
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by  assuming  that 

0-^)-^'^i  ^^^ilC-l+xO^^Q^H^O-^)  (2.79) 


Eq.(2.78)  reduces  to 


On 

^Ao    ,^^  L-^t'  ^^-  y<,^^^^-'       ~^^  ^^.  ,        (2.80) 


Straight  forward  differentiation  shows  that  S^P^  is  minimum 

■       -    Ao 

when 


X^=l  ,  ^.?^\^o. 


(2.81) 


From  (2.79)  and  (2.81)  the  analysis  is  valid  provided 
that  ^'^p  +  "^/O  >/ 1.  which  is  clearly  satisfied  in 
practice. 

From  (2.80)  and  (2.81)  the  minimum  sensitivity  to 
variations  in  Ao  is  derived  as 


SaJ^'z  ^^^   /Ao  '         (2.82) 


49 


The  corresponding  value  of  CJ^P^  is  given  by 


OW 


Ao   - 


AoQp 


(2.83 


LP  realizations  as  circuits  1  and  2  can  be  obtained  by  using 
a  conversion  function  of  the  form  or  its  reciprocals.  The 
admittances  Yl  to  Y8  are  chosen  as 


i^-'(^\    •^2z?,C2.U.S.Cl>-&3    .  \'^r(^A,  V5-0-3,  V/,,>(^-_ 


O  >\V.,  Gh      (2.84) 


in  order  to  obtain  a  conversion  function  of  the  form  V.^--'*- ----fS 


Now  Qpa  and  Cv^pa  are  obtained  as 


O,  .  ^>i  'tOt'^^/c-0^t^"V^,)Mol  'i  -1  ^-(^t^--^/,^)//-,  I 


7'L' 


^Ava  - 


'  r  ^'■-^^^l/^}(if^-'l^,)^o^(^^^l<i<t)(M(^  l(s-i(.i  AO 


(2.85) 


■J     .      ■^  ^7\ 


^h 


(2.86) 


1h 


(2.87) 
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The  sensitivity  of  Qpa  with  respect  to  variations  in  Ao 
can  be  minimized  following  the  approach  used  earlier.  It  is 
found  that  for  minimum  sensitivity [ 31] 

G^-(5^  ^     (^\zQif  (35  .  (2.88) 


The  minimum  value  of  S^^^  can  be  shown  to  be 

Ao 


Sao"_-  4^  /Ao  (2.89) 


and  the  corresponding  value  of  ^^^  is  given  by 


C  ''^    I 


(2.90) 


The  above  sensitivity  analysis  can  be  extended  to 
realizations  using  any  other  type  of  conversion  function 
[31]. 

Equation  (2.69)  shows  that  the  sensitivities  to  passive 
element  variations  are  independent  of  the  selectivity. 
Furthermore,  the  sensitivities  with  respect  to  variations  in 
the  amplifier  gain  are  low.   The  proposed  realizations  are 
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seen  to  have  similar  sensitivity  properties  as  the  low 
sensitivity  realizations  reported  in  [21-27]. 
E.   STABILITY 

It  has  been  shown  elsewhere  [20]  that  some  networks  using 
QIC's  can  be  conditionally  stable  where  a  circuit  can  lock  in 
an  unstable  mode  during  activation  (just  after  switching  on 
the  power  supply) .  In  this  section  the  stability  properties 
of  the  configuration  show  in  Fig.  (2.2(c))  are  examined. 

The  natural  frequencies  of  the  circuit  in  Fig. (2.3(c)) 
are  the  zeros  of  the  characteristic  polynomial  D9s)  as  given 
by  (2.72).  The  differential  open-loop  gain  of  a  frequency 
compensated  Op.  Amp.  ^  in  a  bounded  frequency  range  C  :rLL£-vi:'N\  :. 

where  Ao  and  wc  are  the  d.c.  gain  and  the  cutoff  frequency 
respectively,  and  0  ^  ^o  ^  L'vi.o;?,  .  In  the  frequency  range 
'D<<'sijc  the  amplifier  gains  Al  and  A2  can  be  assumed  to  be 
real.  For  any  second-order  transfer  function  the  coefficients 
of  D(s)  are  seen  to  remain  positive  for  any  attainable  pair 
of  Al  and  A2 .  this  is  due  to  the  absence  of  negative  terms 
in  D(s).  Therefore,  the  zeros  of  D(s)  will  remain  the 
left-half  s-plane  and  low  frequency  unstable  modes  cannot 
arise  during  activation. 
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III.   PROGRAMMABLE  GIC  FILTER 

A.   GENERAL 

Signal  processing  devices  evolved  considerably  over  the 
last  several  years.  The  progress  was  motivated  by  the 
advancement  in  film  and  semiconductor  technologies,  as  well 
as  the  continuous  upgrading  of  systems  specifications  to  take 
advantage  of  the  available  technologies  to  the  limits. 

Linear  filtering  finds  many  applications,  such  as  speech 
processing  (recognition  or  synthesis),  geology, 
instrumentation,  communications,  process  control,  adapting 
balancing,  etc.  There  has  been  much  emphasis  on  performing 
the  filter  function  digitally,  largely  because  of  the  ease  of 
varying  and  optimizing  the  transfer.  However,  and  for  many 
reasons,  such  as  cast  size  signal  processing  complexity,  and 
bandwidth,  it  would  be  desirable  to  perform  the  filter 
function  with  linear  components,  yet  retain  the  flexibility 
of  varying  the  filter  parameters  digitally. 

Recently,  several  advantages  of  combining  linear 
components  (amplifiers  and  capacitors)  and  nonlinear  elements 
(switches)  have  been  demonstrated  using  MOS  switched 
capacitor  techniques  [31,  32].  Here,  we  are  presenting  the 
results  of  realizing  a  continuous  active  device  using  linear 
elements  and  switches  controlled  by  digital  signals  to 
achieve  fully  programmable  filters. 
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Our  research  addressed  two  different  aspects  of 
programmability  namely; 

(1)  Programming  the  filter  topology  using  a  minimal  set  of 
elements  to  obtain  any  type  of  filtering  function  desired, 
e.g.,  LP ,  HP ,  BP ,  N  and  AP . 

(2)  Programming  the  filter's  transfer  function  parameters, 
(pole  resonant  frequency  a^p  and  quality  factor  Qp)  for  a 
chosen  type  of  filtering  function. 

B.   THE  PROPOSED  GIC  PROGRAMMABLE  FILTER 

The  basic  active  network  considered  as  the  heart  of  the 
GIC  programmable  filter  is  the  GIC  structure  [31]  of 
Fig. (3.1),  whose  superior  performance  was  established  in  the 
literature  [10,32].  The  filter  transfer  function  was  derived 
using  loop  analysis  in  Chapter  II. 

Table  (3.1),  illustrates  that  for  any  of  the  LP,  HP,  BP, 
N  and  AP  realizations,  five  resistors,  two  capacitors  and  two 
Op.  Amps.,  are  required,  also,  the  transfer  function  of  each 
realization  is  shown.  The  passive  elements  are  connected  to 
the  different  nodes,  shown  in  fig. (3.2),  for  the  different 
realizations.  A  set  of  MOS  bilateral  switches  controlled  by 
a  digital  binary  word,  are  used  to  interchange  the  elements 
to  achieve  the  different  types  of  filter  realization  shown  in 
Fig. (3.3).  The  truth  table  of  the  switch  control  logics  is 
shown  in  Table  (3.2).  Fig. (3.3(a))  illustrates  the  CMOS 
logical  circuit  for  realizing  this  truth  table.  While  four 
of  the  resistors  are  equal  and  of  value  R  each,  the  fifth 

54 


l-.i. 


3.1   -    The  Generalized  ImniCtance  Converter  (GIC) 
Implementation  Using  OA's 


S€f, 


Filter 

Type 

\ 

^2 

''3 

\ 

''5 

.  ^6 

^7 

^8 

Transfer  Function 

LP 

G 

C 

^^ 

G 

G 

0 

0 

G 

T2  =  2Wp^/D(s) 

HP 

G 

G 

C 

G 

0 

G 

C 

G 
QP 

T^  =  2S^/D(s) 

BP 

G 

G 

C 

G 

0 

G 

G 
QP 

C 

T^  =  2(Wp/Qp)^/D(s) 

N 

G 

G 

c 

G 

G 

0 

C 

G 
QP 

T„  =  (S^+W  ^)/D(s) 
2        n 

AL 

G 

G 

c 

G 

G 

0 

c 

G 

QP 

Tl=(3  -  Qp  s  +wp'')/D(< 

where  T(s)  =  N(s)/D(s)  and  D(s)  =  S^  +  (Wp/Qp)S  +  Wp^ 


TABLE  3.1.  -  The  Elements  Identification  for  Different 
Realizations  of  the  GIG  Filter  . 
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Fig.   3.2  - 


Schematic  Diagram  of  the  Programmable  GIC 
Filter  Showing  the  Controlled  Nodes 
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Fig.  3.3  -      Different  Elements  Realizations  and  the 

Corresponding  Switches  Used  for  Digitally 
Selecting  the  Filtering  Type  . 
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TABLE   3  .2  -The  Truth  Table  of  the  Switches  Logic  Used  to 
Select  the  Filtering  Function 
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Fig.  3.3a  The  CMOS  Logic  Diagram  Used  to  Control  the 
Analog  CMOS  Switches  of  Fig.  and  to 
Realize  the  Truth  Table 
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resistor  is  the  Qp  determining  resistor  and  of  value 

Rq=RQp.   The  two  capacitors  are  equal  and  of  value  C=l/WpR. 

each.   The  two  equal  banks  of  capacitors  are  used  to  control 

ojp.    Each   bank   contains   n   binary   weighted   capacitors 

connected  in  series  through  analog  CMOS  switches  as  shown  in 

Fig. (3.5).   Using  a  digital  binary  word  of  n  bits  to  control 

Wp,  2^  different  values  of  C  will  result  at  the  2  terminals 

of  both  capacitors  banks  that  correspond  to  2^  different 

values  of   p.   Using  a  similar  technique  the  value  of  Rq  can 

be  controlled  through  a  bank  of  m  binary  weighted  resistors 

in  series,  though  analog  CMOS  switches  as  shown  in  Fig. (3.6). 

Using  a  digital  binary  word  of  ra  bits  to  control  Qp,  2"^ 

different  values  of  Rq  can  be  achieved  that  correspond  to  2 

different  values  of  Qp.   Thus,  full  independent  control  of 

the  pole  pair  UOp  and  Qp  are  achieved  by  programming  the 

switches  to  obtain  the  corresponding  C  and  Rp.   It  can  be 

easily  shown  that  with  minor  modifications,  an  additional 

programmable  element  can  be  added  for  the  control  of  the 

notch  frequency. 

C.   THE  REALIZED  GIC  PROGRAMMABLE  FILTER 

A  complete  circuit  diagram  of  the  constructed  GIC  filter 

is  shown  in  Fig.  (3.4).   The  values  of  m  and  n  were  selected 

to  m=n=4.   Thus,  15  different  values  of  y^p  (fp)  and  Qp  were 

obtained  as  it  is  illustrated  at  the  corresponding  Table 

(3.3)  and  (3.4).   the  (designed)  banks  of  the  resistors  for 

the  control  of  Qp  and  the  capacitors  for  the  control  of  oOp 

(fp)    along   with   their   control   switches   are   shown 

correspondingly  in  Fig .(3.7)  and  Fig .(3.8). 
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Fig.  3.4  - 


The  Complete  Circuit  Diagram 
of  the  Programmable  QIC  Filter 


62 


digital   word   controlling  Wp 
I ) 


r 

^0 


So-' 


—    * 


Ci~      s, 


r 


'2 


s:- 


^ ^ 


0 


n 


::>/ 


-Co 


r. 


c         — 


c, 


TV" 


^;- 


-/y        — 


r 


Fig.      3.5    -  The    Two    Capacitor    Banks    Realizations    for    the 

Programming    of    cj 
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Fig.   3»6  -      The  Resistor  Bank  Used  to  Realize  R  Needed 
for  the  Programming  of  Q  . 
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switch 
control 

C 

Fp 

Sa 

Sb   Sc 

Sc3 

nF 

Khz 

0 

0      0 

0 

5.96 

16.693 

0 

0      0 

1 

6.34 

15.709 

0 

0      1 

0 

7.  16 

13.887 

0 

0     1 

1 

7.72 

12.875 

0 

1      0 

0 

8.42 

11.810 

0 

1      0 

1 

9.20 

10.856 

0 

1      1 

0 

10.00 

10.043 

0 

1      1 

1 

11.00 

9.047 

0      0 

0 

15.  10 

6.650 

0      0 

1 

17.60 

5.637 

0      1 

0 

20.80 

4.823 

0      1 

1 

26.00 

3.828 

1      0 

0 

35.50 

2.805 

1      0 

1 

55.00 

1.810 

1      I 

0 

100.00 

.995 

Table  3.3   The  four  -  bit  worcis  that  control 
Fp  anc3  the  corresponding  capacitor 
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switch 
control 

Rq 

Qp 

Sa 

Sb   Sc 

S<3 

K 

0 
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0 

24.0 

15 

0 

0    0 

1 

22.4 

14 

0 

0     1 

0 

20.8 

13 

0 

0     1 

1 

19.2 

12 

0 

1     0 

0 

17.6 

11 

0 

1      0 

1 

16.0 

10 

0 

1     1 

0 

14.6 

9 

0 

1     1 

1 

12.8 

8 

0     0 

0 

11.2 

7 

0     0 

1 

9.6 

6 

0     1 

0 

8.0 

5 

0     1 

1 

6.4 

4 

1     0 

0 

4.8 

3 

1     I 

1 

3.2 

2 

1     1 

0 

1.6 

1 

Table  ? .  4   The  four  -  bit  worcSs  that  control 
Rp  and   Qp 
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IV.   COMPUTER  SIMULATION  OF  THE  PROGRAMMABLE  FILTER 

A.   INTRODUCTION 

In  order  to  observe  the  theoretical  responses  of  the 
different  realizations  as  well  as  to  compare  them  with  the 
experimental  measurements,  two  computer  programs  were  written 
in  Fortran.  those  programs  are  shown  in  Appendices  A  and  B. 
The  first  one  simulates  the  programmable  filter's  responses 
Tl(s),  T2(s)  and  T3(s)  as  functions  of  the  network's 
admittances  Yi,  i=l  .  .  .8,  and  the  "constant"  Op. Amps,  gains 
Al  and  A2.  A  realistic  value  of  A2=A2=10  was  given  to  the 
above  gains  corresponding  to  741  Op.  Amps,  used  in  the 
experiment.  the  second  program  simulates  the  programmable 
filter's  responses  Tl(s),  T2(s)  and  T3(s)  as  functions  of  the 
network's  admittances  Yi,  i=l  ...  8  and  the  frequency 
dependent  Op.  Amps,  gains  Al  and  A2 .  In  this  second  case  a 
single  pole  approximation  value  of  wils  was  assigned  for 

Ai(i=l,  2).   Where,  Wi  is  the  GBWP  of  the  Op.  Amps,  used,  a 

g 
value  of  Wl=W2=2-:t  xlO   was  given  to  the  above  Wi(i=l,  2) 

corresponding  to  741  Op.  Amps,  used  throughout  the  research. 

To  be  able  to  compare  the  computer  simulation  results 

with  the  experimental  ones  obtained,  the  same  values  of 

admittances  were  used  as  input  data.   That  means  that  the 

same  values  of  R,  Rq  and  C  according  to  Chapters  III  and  V 

were  used.   With  the  value  of  R  selected  to  be  1.6k   Table 
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(4.1)  illustrates  the  different  values  of  Rqs  used  in  the 
simulation  to  realize  the  different  values  of  Qp  and  Table 

(4.2)  illustrates  the  different  values  of  capacitors  used  in 
the  simulation  to  realize  the  different  frequencies  Wps. 

The  "DO  CASE  I"  command  of  Fortran  simulated  the  digital 
logic  (including  the  control  switches)  used  to  realize  the 
different  responses  of  the  programmable  filter  that  is  LP, 
HP,  BP,  N  and  AP.  The  frequency's  dx^p)  translations  and  the 
different  QP  values  were  simulated  by  changing  the  values  of 
Rq  and  Cs  at  every  run  of  the  program. 

The  above  programs  also  simulated  the  transfer  function 
of  two  cascaded  GIC  programmable  filters  as  it  will  be 
discussed  later  in  Chapter  VI.  Each  of  the  filters  could 
have  been  at  different  realization  (as  well  as  at  different 
p  and  Qp)  relative  to  the  other  one.  The  25  possible 
combinations  of  transfer  function  realizations  are  shown  in 
Table  4.3. 
B.   SIMULATION  RESPONSE(S) 

1 .   Low  Pass  (LP)  Realization 

Using  the  elements  values  prescribed  in  Table  (3.1) 
yields  the  following  L.P.  transfer  function: 

\2Cb)r (4.1) 


70 


For  R=1.6k« 

Rq 

QP 

1.1 

0.7 

1.6 

1.0 

3.2 

2.0 

6.4 

4.0 

11.2 

7.0 

16.0 

10.0 

32.0 

20.0 

Where  Qp  =  R^ 


Table  4.1  -  The  Resulting  Qp  Values  for  the  Different 
Introduced  Rq  Values  in  the  Computer 
Simulation. 
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For  R-1.6k2 


c 

fp 

nf 

khz 

100.00 

.099 

50.00 

1.99 

26.00 

3.83 

12.5 

7.96 

6.6 

15.10 

3.3 

30.16 

0.85 

117.09 

0.42 

236.67 

fp  = 


2  K-   R.C, 


Table  4.2  -  The  Resulting  tup  Values  for  the  Different 
Introduced  C  Value  (R  constant  1.6k^ 
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fiite; 
•  IP 

'HP 
[BP  , 
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LP 

LP 

LP 

I? 

HP 
-   HP 

HP 

H? 

BP 

ap 

3P 
BP 

a 

N 
H 

il 

AP 
AP 
AP 


TtiVjle.  ^ !)   -   The  25  Possible  Combinations  of 
i  '        Transfer  Function  Realizations 


FILTER  2 

J 

T 

LP 

2 

2 

H£ 

1 

1 

BP 

a 

4 

N 

3 

3 

A? 

3 

3 

HP 

2 

1 

3P 

2 

a 

N 

2 

3 

AP 

Z 

3 

LP 

1 

2 

UP 

1 

^ 

l^ 

1 

J 

A? 

1 

3 

LP 

M 

2 

HP 

4 

1 

IS 

a 

3 

AP 

■  a 

3 

.IF 

3 

2 

BP 

3 

1 

BF 

J 

4 

AP 

3 

3 

IP 

3 

2 

H? 

3 

1 

BP 

3 

H 

N 

3  . 

3 
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This  transfer  function  assumes  the  following  complex 
value 

TzCjOiOr;  -QjQp  (4.2) 

with  magnitude  of 


1  l2LiaJp)\  =  2Qpc  ^0\o<^\2G?^\  ^  ab  (4.3) 


Fig.  4.1  illustrates  the  theoretical  "ideal"  LPF 
magnitude  response  for  fp=3.8Khz  and  for  3  different  values 
of  Qp.   The  simulation  results  match  the  equations  of  (4.2). 

Fig.  4.2  illustrates  the  theoretical  "ideal"  LPF 
magnitude  response  for  Qp=2  and  for  3  different  values  of  (-^-p. 
the  Op.  Amps,  gains  Ai,  (i=l,  2)  are  frequency  depended. 
This  dependance  affects  the  magnitude  of  the  filter  and 
causes  a  frequency  shift  from  the  theoretical  value  of  the 
ideal's  case.  Figs. (4.3)  and  (4.4)  illustrates  the  ideal  vs. 
nonideal  theoretical  LPF  amplitude  responses. 

Data  extracted  from  Figs. (4.3)  and  (4.4)  are 
illustrated  in  Tables  (4.4(a))  and  (4.4(b))  simultaneously. 
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LP.F  AMPLRESPONSE  (P=3.828KH) 


M- 


03 

Q 

H 

Q 
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i 

a 


LEGEND 
=  .707 


I        I      r    I    I  I  I  I 


I       I      t    I    I   M  f 


icr 


I I r I I 1 r r| 


FREQUENCY  HZ 


Fig.  4.1  -  "Ideal"  LPF  Amplitude  Response 
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LRF  AMPLRESPONSE  (Q=2) 


m 

Q 

U 

Q 


a. 


LEGEND 
F  =  IKH 
F  =  4KH 


l(f 


I  I  I        I       I     I     I    r 


FREQ 


Fig.  4.2  -  "Ideal"  LPF  Amplitude  Response 
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LP.F  AMPLRESPONSE  (Q=2) 


FREQ 


ideal 
nonideal 


Fig.  4.3  Ideal  .vs.  nonideal  L.P.F.  amplitude  response  for  Q-2 
and  frequencies  (1.99K,  7.96K,  30.16K  and  117.lK)h> 
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LRF  AMPHTUDE  RESPONSE 


10* 
FREQUENCY 


1. 

Q   =    0.7 

2. 

Q   =    2.0 

3. 

Q   =    4.0 

4. 

Q   =    7.0 

Fig.  4.4   Ideal  .vs.  nonideal  L.P.F.  amplitude  responce  for 
frequencies  {7.96K,  30.1K)Ktand  variety  of  Qs . 
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Data  From  Fig.  (4.4)  LPF  Characteristic 


Ideal                              Nonideal 

QP 

|T2(s)| 
at  pick  value 

Rq 

QP 

|T2(s)| 

at  pick  value 

1 

1 

.07 
2.0 
4.0 
7.0 

3.0 
12.0 
18.1 
22.9 

1.1 

3.2 

6.4 

17.2 

6.7 
2.0 
4.0 
7.0 

3.0 
12.0 
18.1 
22.7 

The  pick  values  obtained 
at  fpick  =  7.9k 

The  pick  values  obtained 
at  fpick  -  7.9k 

1 

(b) 

Ideal 

1 

Nonideal 

Fp 

JAmp 

.  Response 

i    ^^ 

fp 

Amp 

Response! 

khz 

. 

dB 

' 

khz 

i 

dB 

1990 

.45 

12. 

04 

50 

1990. 

45 

12. 

72 

7165 

.60 

12. 

24 

12. 

5 

8359. 

87 

12. 

113 

28662 

.42 

12. 

31 

6. 

6 

25875. 

.80 

12. 

04 

109745 

.50 

12. 

31 

0, 

085 

90565. 

00 

11. 

08 

Table  4.4  -  Data  from  Figs.  (4.3)  and  (4.4)  indicating 
the  Affect  of  Frequency  Dependency  of 
Ai,  i=l,  2 
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2.   High  Pass  (HP)  Realization 

Using  the  elements  value  shown  in  Table  (3.1)  yields 
the  following  HP  transfer  function: 


lics^)^      ^^^  (4.4) 


takes   the   following  complex  value   at  wp. 

TiCjOip)  -"^jQ?  (4.5) 

which  magnitude   is 

U^Uoop)  l-2Qpc  2Dlo«3.C:ja)p^  ;  ab  (4.6) 


Fig.  (4.5)  illustrates  the  theoretical  "ideal"  HPF 
amplitude  response  for  f p=3 . 8KHZ  and  for  3  different  values 
of  Qp.  The  computer  simulation  results  match  that  of  the 
(4.6)  relation. 

Fig.  (4.6)  illustrates  the  theoretical  "ideal"  HPF 
amplitude  response  for  Q=2  and  for  3  different  frequencies. 

Fig.  (4.7)  illustrates  how  this  dependency  effects 
the  amplitude  of  the  HP  filter.   (The  amplitude  decreases  as 
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H.P.F  AMPL  RESP0NSE(P=3.828KH) 


o 
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Q 
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Q 
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o 

2 
I 


LECSND 
^  /  Q  =.707 

3/^  =  7. 


■n 
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I     I    M   I  I 
Irf  l& 

FREQUENCY  HZ 


I  I  I  1 1 


Itf 


»     I    I   I  1  r 


icT 


Fiq.  4.5  -  "Ideal"  HPF  Amplitude  Response 


for  a  Variation  of  Qs . 
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H.P.F  AMPLITUDE  RESP0NSE(Q=2) 


n 

Q 
Q 
H  -■ 

a 

flu 

a 


«- 


V 


I  I  I  1 1 1 1 

FREQUENCY  HZ 


LEGEND 
=  1KH 

=  4KH 
F  =  9KH 


I  '  '  I  f 


itf 


Itf 


Fig.  4.6  -  "Ideal"  HPF  Amplitude  Response. 
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H  J.F  AMPUTUDE  RESP0NSE(Q=2) 


LEGEND 
=  4 


H= 


FREQ 


Fig.  4.7  -Nonideal  HPF  Amplitude  Response 


83 


the  frequency  increases)  while  Figs.  (4.8)  and  (4.9) 
illustrate  the  frequency  shift  of  the  amplitude  response  due 
to  that  dependency. 

From  the  date  of  Fig.   (4.8),   Table   (4.5)   were 
constructed. 

3.   Band  Pass  (BP)  Realization 

Using  element  values  from  Table  3.1  the  following  BP 
transfer  function  is  achieved: 


IWM-  ^^Q?^^ (4.7) 


takes  the  following  value  atujp, 


T'lU^i^P^zQL  (4.8) 


which  has  constant  magnitude  of  6dB 


lTiijvV>p)U2^  lOlo^Z  =  ^<ib  (4.9) 
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H.P.F  AMPUTUDE  RESP0NSE(Q=2) 


«>-i 


FREQ 


Ideal 
nonideal 


Fig.    4.8    -    "Ideal"    vs . Nonideal   HPF   Amplitude 
Response . 


Ideal 


Nonideal 


fp        Amp  Response 

nf 

fp 

Amp  Response 

Hz 

dB 

Hz 

Db 

4,378.98 

12.28 

26 

4,378.98 

12.148 

8,359.87 

12.30 

12.5 

8,359.87 

12.11 

15,923.57 

12.31 

6.6 

15,923.57 

12.02 

31,847.13 

12.31 

3.3 

31,847.73 

11.41 

123,407.00 

12.31 

0.85 

107,484.00 

9.21    i 

252,781.05 

12.31 

0.43 

187,101.80 

9.21    i 

242,781.05 

1 

12.32 

0.43 

187,101.80 

6.69 

i 

Table  4.5  -  Data  Illustrating  the  Ideal  vs.  the  Nonideal 
Responses  of  the  HPF 
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H.P.F  AMPL  RESP0NSE(P=8K/16K) 


EQUENCY  HZ 


1.  Q   -    0.7 

2.  Q   =    2.0 

3.  Q   =    4.0 

4.  Q   =    7.0  ^         ,  . 

Fig.  4.9  -  "Ideal"  vs.  Nonideal  HPF  Amplitude 
Response . 
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Fig.  (4.10)  illustrates  the  theoretical  BPF  magnitude 
response  for  fp=3.83  KHZ  for  different  values  of  Qp.  These 
agree  with  the  (4.9)  equation  since  as  it  is  indicated  by  the 
simulation  plot  the  amplitude  is  constant  and  independent  of 

Qp. 

Fig.  (4.11)  illustrates  the  above  concept  but  at 
fp=15.1KHZ. 

Fig.  (4.12)  and  (4.13)  illustrate  the  theoretical 
"ideal"  BPF  amplitude  response  for  Qp=2  and  for  different 
values  of  p.  As  it  is  indicated  from  Fig.  (4.13)  the 
amplitude  remain  constant  even  at  very  high  frequencies  (10 
HZ).  But  with  Ai,  (i=l,  2)  depending  on  frequency  the 
amplitude  decreases  as  the  frequency  increases.  This  is 
indicated  in  Figs.  (4.14)  and  (4.15)  which  describe  the  BPF 
amplitude  response  plots  for  Q=4  and  Q=l  respectively  and  for 
different  frequencies.  The  frequency  dependence  of  Ai,  (i=l, 
2)  creates  a  frequency  shift  from  the  ideal  theoretical  value 
which  is  indicated  in  Figs.  (4.16)  and  (4.17).  Table  (4.6) 
illustrates  the  data  extracted  from  figs.  (4.16)  and  (4.17). 
4.   Notch  (N)  Realization 

Using  the  admittances  value  of  Table  (3.1),  the 
following  Notch  transfer  function  can  be  achieved: 


(4.10) 
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B.P.F  AMPL  RESP0NSE(P=3.828KH) 


5*l(f 


FREQUENCY  HZ 


Fig„  4.10  -  "Ideal"  BPF  Amplitude  Response. 
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B.P.F  AMPLITUDE  RESP0NSE(P=16K) 
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Fig.  4.11  -  "Ideal"  BPF  Amplitude  Response. 
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B.P.F  AMPLITUDE  RESP0NSE(Q=4) 
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Fig.  4.12  -  "Ideal"  BPF  Amplitude  Response 
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B.P.F  AMPLITUDE  RESP0NSE(Q=4) 
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Fig.  4.13  -  "Ideal"  BPF  Amplitude  Response 
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B.P.F  AMPLITUDE  RESP0NSE(Q=4) 
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Fig.  4.14  -  Nonideal  PBF  Amplitude  Response 
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B.P.F  AMPLITUDE  RESP0NSE(Q=1) 


Fig.  4.15  -  Nonideal  BPF  Amplitude  Response. 
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B.P.F  AMPLITUDE  RESP0NSE(Q=4) 
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Fig.    4.16    -    Ideal    vs.    Nonideal    BPF   Amplitude 
Response . 


95 


B.P.F  AMPLITUDE  RESP0NSE(P=16K) 


Ideal 
nonideal 


Fig.    4.17    -    Ideal    vs.    Nonideal    BPF   Amplitude 
Response . 
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Ideal 

C 

Nonideal 

1 

fp         Amp.  Response 
khz            dB 

nf 

fp       Amp.  Response 
khz          dB 

7,961.78        6.02    j      12.5 
30,254.78        6.02          3.3 
115,445.80        5.97          0.85 

1 

7,969.738      5.958    , 
28,662.42       5.75     ! 
95,541.38       4.92     j 

• 

— . . 

and 


Ideal 

Rq 
k 

Nonideal 

QP 

Max. 
Amp.  Response 
dB 

QP 

Max. 
Amp.  Response 
dB 

4 
8 
20 

6.01 
6.01 
5.95 

6.4 
12.8 
32.0 

4 
8 
20 

5.902 
5.728 
5.528 

fp^  =  15, 127. 33Hz 

fPl 

=  14, 739. 29Hz 

Table  4.6  -  Data  Indicating  the  Ideal  vs.  the  Nonideal 
Responses  of  the  BPF  Realization 


97 


N.F  AMPUTUDE  RESP0NSE(Q=2.) 


2*lCf 


FREQUENCY  HZ 


Fig.  4.18  -  Ideal  Notch  Amplitude  Response 
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NOTCH  AMPLRESP0NSE(P=3.83KH) 


FREQUENCY  HZ 


Fig.  4.19  -  Ideal  Notch  Amplitude  Response 
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N.F  AMPUTUDE  RESP0NSE(Q=2.) 


md(f 


FREQUENCY  HZ 


Ideal 
nonideal 


Fig.    4.20    -    Ideal   vs.    Nonideal    Amplitude    Response 
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where  n  is  the  Notch  frequency. 

At  P  the  transfer  function  takes  the  value 

TaCjo.,) -.  J^^Ll^^  (4.11) 

Fig. (4.18)  illustrates  the  "ideal"  theoretical  Notch 
filter  amplitude  response  for  a  variety  of  frequencies  and 
constant  Qp(Q=2),  while  Fig. (4.19)  for  a  variety  of  Qps  and 
for  constant  p  (fp=3 .83KHZ) . 

Fig.  (4.20)  illustrates  the  effect  of  the  frequency 
dependency  of  Ai,  (i=l,  2)  on  the  response. 
5.   All  Pass  (AP)  Realization 

As  proposed  in  Table  (3.1)  using  the  same  admittances 
values  as  Notch.  An  All  Pass  transfer  function  can  be 
derived  as 


\\LS'):=       ""    "^P 


.^.;^s>^p^ 


S^-v^ets^^V 


(4.12) 


9? 
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which  takes  the  following  values 


o 


->  O 


(4.13) 


The  above  agree  with  the  computer  simulation  results 
of  Figs.  (4.21)  and  (4.22). 
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A.P  AMPLITUDE  RESP0NSE(Q=2) 
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Fig.  4.21  -  Ideal  vs.  Nonideal  APF  Amplitude 
Response. 
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A.P  PHASE  RESP0NSE(q=2/P=8K) 
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Fig.  4.22  -  Ideal  vs.  Nonideal  Phase  Response. 
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V.  REALIZATION  OF  PROGRAMMABLE  GIC  FILTER 

A.   EXPERIMENTAL  RESULTS 

After  the  circuit  of  Figs.  (3.10),  (3.9),  (3.8),  (3.7), 
was  constructed  a  variety  of  measurements  were  taken  in  order 
to  study  the  response  of  the  network  to  the  different  inputs 
(control  bitwords).  To  observe  the  affect  of  the  control 
switches  which  introduce  a  resistance  of  80  each  at  CLOSED 
position,  two  values  of  Rs,  Rq,  and  Cs  were  used  with  one 
decade  difference  in  magnitude.  That  means  that  R  was  given 
the  values  of  1.6K  (as  discussed  in  Chapter  III)  and  16K,  the 
four  resistors  that  consisted  the  Rq  bank  were  of  values 
(1.6K,  3.2K,  6.4K,  12. 8K)  in  the  first  case  and  (16K,  32K, 
64K  and  128K)  in  the  second  one,  and  that  the  capacitor 
bank's  capacitors  were  chosen  of  values  (lOOnF,  50nF,  12, 
5nF,  llnF)  and  (lOnF,  5nF,  1.2nF,  O.lnF)  accordingly,  to  be 
able  to  keep  the  range  of  frequencies  as  much  the  same  as 
possible  for  both  cases. 

1.   Low  Pass  Filter 

With  the  topology-control  bit  word  000  the  network 
realized  a  LPF  response.  Fig.  5.1(a)  illustrates  the  LP 
response  for  R=1.6K'ito  a  variety  of  frequencies  for  Q=5 
while  5.1(b)  illustrates  the  same  but  for  Q=2.  It  can  be 
observed  from  5.1(a)   (Q=5)   that  the  amplitude  response 
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(b)  LPF  Response  for  R=1.6k  (Q=2,  Rq=3.2k; 
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decreases  while  the  frequency  increases  as  it  was  expected 
from  computer  simulation.  This  is  occurred  up  to  the 
frequency  of  6KHZ;  then  it  started  increasing  with  the 
frequency,  while  at  Fig.  (5.1(b))  (Q=2)  it  remained  constant. 
Fig.  (5.2)  illustrates  (for  R=16KHZ)  that  the  observation  in 
Fig  (5.1)  is  not  any  more  the  case  (for  that  frequency  range) 
and  the  network  responses  the  same  as  in  computer  simulation 
while  in  Fig.  (5.3)  (Q=5)  the  above  can  be  noticed  again, 
this  is  due  to  the  interference  of  the  control  switches  as  it 
was  discussed  previously.  Figs.  (5.4)  and  (5.5)  illustrate  a 
variation  of  Q  values  for  f=9KHZ  for  the  two  cases  (R=1.6k 
and  R=16K  )  respectively.  A  difference  in  magnitude  can  be 
observed  due  to  the  interference  of  the  control  switches. 
Figs.  (5.6)  and  (5.7)  illustrate  the  same  but  for  F=12.8KHZ. 
Figs.  (5.8)  and  (5.9)  illustrate  the  phase  and  amplitude 
response  of  the  LPf. 

2.   High  Pass  (HP)  Realization 

With  the  topology-control  bitword  001,  the  network 
realized  a  high  pass  filter.  Figs.  (5.10)  and  (5.11) 
illustrate  the  HPF  amplitude  response  for  a  variety  of 
frequencies  and  (Q=2).  It  can  be  observed  (as  in  LPF 
realization)  that  for  this  frequency  range  and  for  R=1.6Rt 
the  amplitude  starts  increasing  as  the  frequency  increases, 
decreases  then  again  while  for  R=16K<£  this  does  not  occur, 
the  same  observation  submerges  comparing  Fig. (5.14)  and 
(5.15).   Figs.  (5.12)  and  (5.13)  are  the  plots  obtained  for 
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f=6.65KHZ  and  a  variety  of  Qs  for  R=1.6K^  and  R=16K9  , 
respectively.  a  difference  of  about  12dB  in  amplitude 
appears.  Finally,  figs  (5.16),  (5.17),  (5.18),  (5.19) 
illustrate  the  phase  and  amplitude  responses  of  the  HP 
realization. 

3.  Band  Pass  (BP)  Realization 

The  010  topology-control  bitword  realizes  the  Band 
Pass  Filter.  Figs.  (5.20)  and  (5.21)  illustrates  the 
amplitude  response  for  set  of  frequencies  and  Q=10.  At  both 
cases  the  amplitude  decreases  while  the  frequency  increases 
until  the  value  of  9KHZ.  Then  starts  increasing  again.  It 
can  be  also  observed  that  for  R=l . 6K  the  frequency  deviation 
from  the  theoretical  fp  is  larger. 

Figs.  (5.22)  and  (5.23)  illustrate  the  amplitude 
response  again,  but  for  Q=l.  This  time  the  deviations  from 
the  theoretical  response  (both  of  amplitude  fluctuations  and 
frequency  shift)  are  less. 

Figs.  (5.24),  (5.25),  (5.26),  and  (5.27)  are  also 
plots  of  the  amplitude  response  but  for  variation  of  Q.  A 
difference  of  approximately  IdB  appears  for  the  lowest 
frequency  (3.8KHZ)  and  of  0 . 5dB  for  the  higher  ( 9KHZ ) , 
Finally,  Figs.  ( 5 . 28 ) - ( 5 . 32 )  illustrate  the  phase  and 
amplitude  of  the  BPF  response. 

4.  Notch  (N)  Filter  Realization 

With  the  topology-control  bitword  Oil,  a  Notch  filter 
realization  can  be   achieved.    Figs.   (5.33)   and   (5.34) 

122 


J<  F  R     F  C 
X  F  R     F  C 


H     M  K  R  ;  2  .  7  d  B 

H  =     0  o     CENTER 


2    dB-'DI'.; 
5  0  -^  .^  D  I  ^V 


\1 


i-ji«- 


1— -.-:}..,_.,^^.^...,-F...-.._. 

•2. 


,..,-v-,__.,^r 


-*7 


LI 


•x     GHz 

n  K  R  ^         1     4  0  Q     H : 


2  5     !•:  H  z 
B  I'J  ^      2  0  0     H : 


CA) 


:^^  F  R 

X  F  R 

FCl 

F  C  1 

■  \\     M  K  R  ■■ 
'\\  ■■      0  o 

C  E  H  T  E  R 

?dB 

■7» 

c 

d  B  ■• 

0  O  .■■ 

D  1  '...' 
D  I  U 

,-■."-, 

''" 

s 

,.-1^ ^ 

■--•■■■->- 

-2. 

1 

t 
1 

•\ 

-.-.. 

-     .- 

,...,....-., 

/  ■.-•'"".   i' 

"''*'s.*"  "l .- 

i' 
1 

""■'--- 

- — — 

---.— 

---■— ---^ 

/ 

/ 

1 

/ 



■  -.     0     Hz 

M  K  R  =         3     3  0  0     H 


2  5     K  H  z     .-• 
B  H  :     2  0  0     Hz 


1 .  Phase  '        • 

2.  Amplitude 

Fig.  5.16  -  HPF  Amplitude/Phase  Response 
(Q=l) 


123 


XFR 
XFR 

FCl 
FCl 

r  H  \1  K  R  : 

1  . 
CENTER 

4dB 

2 

dB.- 

0  O  ./ 

DI  '.; 

DI '.; 

r- 

'-».. 



■■ 

-_4. 

____ 

'Z, 

-•^""' 

~— 

- 

— — - 

1 

/ 



1 

GHz 
M  K  R  ■■         h 


0  0   H2 


2  5  K  H  T 
B IJ  =  2  0  0  H 


Y.  F  R 
.■^  F  R 

FCl 
FCl 

■  W     f1  f(  R  ■■ 
-\\-      0  o 

2 
CENTER 

3dB 

c 

dB- 
0  o  ,• 

DI  '.; 
D  I  '.y 

V 

'-I 

>. 

2. 

.M-.-. 

.■•.  ------ 

---•.  .•-•- . , 

-— — .- 

..__ 

. 

/ 

\ 

••.  0  H  z 

M  K  R  ^    2  4  0  0  H 


2  5  f- :  H  z 
B  N  •■     2  0  0  H 


1.  Phase 

2.  Amplitude 


Fig.  5.17  - 


HBF  Amplitude/Phase  Response 
(Q=0.5) 


124 


XFR 
XFR 

FCl 
FCl 

r  N    M  f- :  R  : 

'W--      0  o 

3  .  I 
C  E  H  T  E  R 

5dB 

dB.' 

0  o  /• 

DI  V 
D  I'^ 

^-'-- 

V-x. 

r''- 
,/ 

"■•-X... 

I 

I 

/ 

y 

■-v-.- 

2. 

r.     ...    , 

> 

^ 

1 

— - 

— ^ 

_._i_ 

. 

1 

1 

••    0     Hz 
M  K R  ■        4     8  0  0     Hz 


2  5     K  H  2 
B  W  ^      2  0  0     H 


1  -1 
i'x 

FR 
FR 

FCl 
FCl 

■  W       M  \\  R   : 

■  N  =     0  o 

C  E  H  T  1  R 

bdB 

d  B .  •■ 
0  o ,.  ■• 

D  I  '.' 
D  I  y 

y     {V 

• 

H 

's.  ^ 

'  'rt 

-2. 

•,-:.\r- 

--,....- 

>  — .i,,-- 

-■■'"'•-" 

•  t*"'"ii\ 

-■•.....-....- 

■-''..i-"-i,. 

j 

»_ 

.- 

d. 

1 



■■■     0     Hz 

M K  R  =         1     2  0  0     Hz 


2  5     Y:  H  z     /' 
E  N  =      2  0  0     Hz 


1.  Phase 

2.  Amplitude 


Fig.    5.18    - 


HPF    Phase/Amplitude    Response 
(Q=2) 


125 


X  F  R     F  C 
X  F  R     F  C 


H     MKR^  3.4dB       2     dB.-DI'.; 

N=     0O     CEhUER  5130/DIV 


^ 


^ 


."^"V 


'=TT- 


2. 


:^»'-;.— j^ 


p;«i  iy,i*>^i^T'  ■( 


n 


^.     0Hz 

M  K  R  •         3     6  0  0     H 


B  U  ■■      2  0  0     Hz 


1.  Phase 

2.  Amplitude 


Fig.    5.19    -    HPF   Amplitude    Response     (Q=2) 


126 


'     i     i 

'     '     t 

1    : 

.  .       \ 

■ 

i   i   i  i 

1  '  1 

M  '  ' 

Mi; 

'    !       1 

.    ■ 

'    ,    i    , 

;     ;     ' 

2dbldi>y      . 

::2y<b2!MiNCM 

!   '  i 

I'll 

i  1  M 

1 

i    M 

i    1    ;    1 

■    1    1 

i  1  M 

i 

;    1 

i    1    i    ' 

!     '     ' 

1     1     : 

!              ' 

1    1    1 

1     1     j 

,     t     f     1 

1    !    1    ^ 

' 

'     '     ; 

MM 

i     ; 

i            ;      i 

:        !        ; 

1    1 

1     1     1 

\    \    '    ' 

1    .  ' 

!       - 

MM' 

i     1     ! 

1    1  ' : 

Mi' 

jTx       i 

MM 

i     1     1     i 

'    : 

'     i     , 

j'"' 

•^O!    ' 

' 

' 

/ 

X 

'  ^1  1        •  ■' 

■       ,       , 

..  ir._^ 

y_  _.„v. 

•i  ' 

1    ; 

,      ■ 

i    !    ■    : 

i? 

-\ 

JU 

1    : 

' 

■    ' 

■    1 

! 

/^ 

■\ 

it 

1    i 

'      :      .      . 

/^f 

_  Z j 

\,    ' 

'  ■  < 

\ 

II 

I    i 

' 

/ 

1  \         / 

■ 

>^ 

~ 

, 

U^A 

Y 

k 

\^ 

^^    / 

; V- 



p— - 

—■ — 

/ 

1 

— yv 

^        / 

—  -  Y 

I      A       7 

1       /     N 

\     /           ' 

/ 

\ 

-1  i 
/ 

1 

/— 

\ 

1    1 

-■-\---l 



1      / ' 

\— — - 

-  /  — 

/— - 

_  1  /-- — 

\ 

k    -.  

/— 

./ 

— \ 

:      ;      1 

Y  1 

1         : 

'    1 

\ 

L]  I 

-  • 

-  r- 

-. . 

^::::i:-\ 

/ 

it' 

1 

/ 

-  — 

„.  \ 

|L  /  ^ 

1          V_ 
/       "    V 

[ 

/— — 

I  Vr      v/r           iv 

L  l 



-.;.::■_ 

tMV  i 

'   '    ■     \ 

~    A '^ 

^— 

— . 

F^l    i    i 

I  ' 

.  T           1  '          1            1           *,  '         1 

.. 



- 

1     - 

f    =    0.99khz 

2    - 

f    =    1.99khz 

3    - 

f    =    3.8khz 

f  =  6.65khz 
f  =  lO.Skhz 
f  =  IS.lkhz 


Fig.  5.20  -  BPF  Amplitude  Response 
(R=1.6k) ,   Q=10 


127 


'^\ — 1 — i — " " 

-  ,       -r~  " 

:     1 

■ 

;    1 

' 

■  ' 

'      ' 

'" 

■ 

1       .       .       : 

'  1  i  1 

I    !    1 

1    !    : 

' 

'    i    '    ' 

: 

idt/ai 

V:     I     ■     . 

! 

1  i  !  ' 

'  '  M 

;    1    , 

i 

!    '    '    ' 

'    ;    !    ' 

:  1  i  ' 

i    i    :    : 

2^rta/4.y     ~^ 

i 

' 

— 

^. 

.       ;      1 

1 

'         :         ■ 

1       '       '      1 

'  1 

'    !    : 

i     1     ■ 

1 

:      1      , 

.    1    , 

-— 

As 

■    '    i    ■ 

JT7(S 

;     ■ 

12;' 

! 

'     ;     ■     ■ 

\i 

», 

11 

/    \ 

' 

' 

i 

U-      _ 

1    /     \ 

-- 

~::\ 

1 

1    1 

y    V 

'    '         1 

f      r 

1 

j 

r      \ 

• 

j 

1 

1  /             \ 

'■ 

1 

1 

T:\f 

=    :fJ 

■ 

-_- 

.^_ 

J 

l^ 

■  • 

.„    T 

-____.^ 

— 

-' 

/ 

1 

— 

■^-n- 

I 

iV-  l:-A 

I 

/ 



1 

,     1     , 

— 

_L_ 

■    -^ 

::z:i: 

- 

\ 

\       ^ 

„.. 

j 

— 

r""  / 

— ■^—- 

-  V-r- 



— 

-    — ^— 

1 

- 1- 

'   \     i 

\ — / 

. — 

j                       / 

\       j 

— 

L_ 

~~\  I 

.— 

y  /  - 

— _-^_., 

^V 

-.-- 

i- 

-■F^ 

1 

- 

—    f 

Y— 

-  -  — —  ■ 

hi               , 

-- 

—  - 

-S 

"^h 

_: 

__ \ 

\" 

1 

1           f 

/         V 

- 

^.^j. 



■ 

1 

.J_:-' — 

,    ; 

:     1 

1.  f   =   0.99K 

2.  f   =    1.99K 

3.  f   =    3.8K 


4.  f   =    7.96K 

5.  f   =11.2    K 

6.  f   =15.1   K 


Fig.     5.21    - 


BPF  Amplitude  Response 
(R=1.6k) ,  Q=10 


128 


'      i 

i    ■    , 

■ 

' 

' 

,      ,      ■      1 

1 

1    i 

i     : 

'   !   i   1 

1    !    , 

'    !    i    ! 

'    i    I    ! 

i    J    !   ' 

'    1    :    ■ 

III; 

[    i    ' 

M    :    I 

i   ■   i   ■ 

i    '    : 

1          1          :         1 

: 

! 

,          •■ 

■  i 

'  i     ' 

1   i 

1          '         * 

1  '  '  ' 

i 

•    •    i 

'   1   ' 

1 

' 

1     i     ' 

1    ■    • 

\   \   ',   '■ 

1 

1     1     ' 

•  i  ,  1 

'     1     1 

1    ;        1 

1     ■ 

1    !    i 

;      ; 

1 

1      I 

'    1    : 

. 

'  \  ■ 

:       1       1 

1     ' 

:      .      i 

!       i        :        , 

:        :    1 

'          I 

, 

i      : 

1     ■ 

' 

I 

a'^ 

1 

1 

3 



^       1 

:  :  ■      4  - 

— -— .^^^ 

l\  p 

2 

J-^ 

, 

"~ 

yC 

__js^ 

_____ 

Y^: 

y\ 

r  -^^^  ,-     !            1 

r":i„        L 

t 

\     ^"^ 

/ 

V- 

\/    5 

Sr-— > 

^ 

/_ 

"V,^    ! 

1 

Y 

A~ 

— . 

^^s^ 

1 

A " 

A — 

'    '    . 

^s,,^           1 

""K^ 

' 

/  \ 

\      : 

/ 

. 

1 

/  V- 

f  A  / 

\     y 

1 

_        .        1^ 

j/\ 

1     : 

1 

/ 

~       ~            I 
_  .    I  _ 

Tp-  )f^ 



k^A 

U_-/^X1 



1 - 

lA 

i\ 

.                    I               .J 

-  .i__,  .  __ 

^  i  \ 

/ 

^^ 

k    - 

f  1  A  J 

/ 

y^ 

/  "    > 

/ 

~| 

^ 

^ 

—     jf  -  - 

— y 

_l 

^v 

— 

— .  — 

*s^     -.  -4     - 

«d 

.-- 

\--T 

/ — 

1 N 

X 

>^ 

1.  0.99khz 

2.  1.99khz 

3.  3.8  khz 


4.  6.65khz 

5.  12.8khz 

6.  IS.lkhz 


Fig.  5.22  -  BPF  Amplitude  Response 
(R=1.6k) ,  Q=7 


129 


Mil 

III! 

:   !   1   I 

1   i   1   ' 

'■        : 

i     ;     '     ' 

— T ! ! ' — ' 

!    !    '     ■ 

1         ; 

~ . 

1   ,      1 

i     1     ' 

1      ' 

'    1    i    ' 

ill! 

M   ' 

1   ,   :   . 

I'll 

i     '     1     ' 

M  a>^v 

! 

i   !   M 

M   ■   1 

;    ;    1 

1 

1 

'  I  !  i 

1   !   1   j 

Mi 

Mm 

!     '     ■ 

'     ! 

:      1      1 

1 

:    ! 

1 

ill! 

i     •     ' 

;    1    ' 

i      '      1 

1 

,    j 

j  '  ■  1 

' 

MM 

i     :  ; 

1      1      ■ 

1  '  '  I 

:       '       • 

;       '       ■ 

Mm 

■ 

!        1    ' 

1             i             i 

'    1     ' 

!    '■    •■ 

\  I  '<  \ 

■ 

Ml! 

•       '       1 

±^\ 

'     ;     i     ; 

Mi 

!    ■    ■ 

1 

i    i    ' 

1 

' 

'   !   1 

i      1      '      ' 

-^  l-\  / 

\^ 

"T^ 

i            ,      1 

!      i      i      :    ~ 

5 

■^ — , — 

^ 

--^^— _ 

,       .      —  — 

-^- 

11 

~\ / 

^,»-v-/_ — ' 

^^"^A^ 

^ 

h><^ 

■^         - 

-r 

7<1 

1   ^"""s^^ 

il  v 

—  V— 

"V — V 

r       ^ 

^ 

——— 

-VS. — 

______ 

—  t  — *\ 

_     l_^ 

r  V 

/—— 

/\" 

^. 

\^ 

/l — \l         \ 

f — ~~ 

y    s 

s^ 

k. 

_4_: 

:^ 

1  ^A —  / 

\.~~  / 

— ^— — 

y\ — 

_.: :_ 

>• 

V 

. 

-K 

.._. 

_ '  ^ 

1^1 

-\r 

\/~ 

_ — y 

""' 

.\L. 

/K~ 

~  / — 

.... 

\v 

>\^ 

T     /      U 

■ 

>v 

^/v, 

/ 

/— 

— - — 

N 

— -- 

—     M^'^v^  r 

'. ' ' 

\r— 

. — ___ 

:-;:---:.  j. 

/    /    Y 

V 

1 

>-  fPP 

\/ 

\— 

X. 

.. . 

1 

\ 

___. 

'"*'"""'''^'"~->~-»              i 

1/ 

1 

— .  1 



1.  0.99khz 

2.  1.99khz 

3.  3.83khz 


4.  7.96khz 

5.  12.8khz 

6.  15.  khz 


Fig.  5.23  -  BPF  Amplitude  Response 
(R=1.6k) ,  Q=l 


130 


1 

i 

1 

1 

1     i          ', 

; 

'     1     . 

1  1  1  i 

, 

Ml; 

1    i       ' 

1     ;     ;     , 

i    i  ' 

'     1     i 

i  M  i 

1 

I  1  i  1 

1    '    ' 

i  '  '  : 

'   !       ' 

'   i   1   ! 

1     1     ; 

1      !  i 

i     ■     i     . 

1 

1  !  '  ' 

'       !       ! 

1    i    ! 

I   . 
1 

Ill 

i       '       , 

1 

1       t       i       ] 

, 

1  '  ' 

1     '     : 

-      ,         i  .    I« 

v;     ' 

1 

'    i    i 

,        :        '        ; 

!     1     . 

5dbM' 

III 

i     !     ' 

' 

!       1 

z^ih^yjcticr: 

i    !    i    : 

; 

1  ] 

:     '     1     ' 

:          1          '          ! 

1  '  1 

1 

'     ■     ' 

:       ■ 

M  :  : 

1 

'     ■     1 

.     ! 

1 

;    ; 

'11 

I 

1     i 

•       1 

■     ' 

: 

-.,—  -, 

i--     3-- 



■     ' 

/^ 

/       > 

/ 

\ 

/ 

\ 



f     . 

\ 

1 

\ 

/ 

\ 



/ 

\ 

-^  / 

'^ 

\ 

k    ■     , 

— -■ 

. , 

Z 

— T 

i          1 

r^ 

\    '     ' 

— —  ■ 

1 

I     \ 

\ 

1 

/       \ 

\ 

^^L 

t-L:^ 

v 

! 



x 

!_,_, — 

__.  ,  -^ 

4-'/A- 

\      '     '     '■ 

i       ,    > 

V 

_ _._._ 

\    ' 

\ 

X 

-\ 

\ 

>V 

1 

\ 

— . 

^  .    .       -  j 

\r-V 

,    • 

/ 

\\ 

^^>>>^ 

r     1  / 

\          \         1 

>N 

1 

K— :-h^ 

-T— : 

-— •  ■- 

-— "H— - 



-  — 

...  „.._.  .„ 

-  :    :  t 

Q 
Q 
Q 


10.0 
4.0 
1.0 


Fig.  5.24  -  BPI  Amplitude  Response 
(R=1.6k) ,  f=3.83khz 


131 


:                 1 

, 

,    ■    ; 

1 

i    !    M 

Mi 

' 

1          , 

Ml 

:     .     1 

1      •      , 

;     .     , 

!   1   !   1 

:     i     i 

MM 

1     1     . 

1    1    1 

!    M 

1     .     ,     ^ 

:      f      :      ■ 

Mil 

i     ' 

1   M 

i     '     ' 

'     : 

I'll 

1   j 

,     ■     , 

, 

1   '   1   ' 

1     ! 

1         ;     ; 

,     ;     1     1 

i      ;      • 

i     j 

'     1     . 

1      i      : 

■      ;      !      1 

< 

1 

' 

1     ' 

1 

■      ,     ; 

M       :      ' 

j 

:    ] 

Mil 

i      ,      • 

■      ' 

1      1 

:     . 

I    , 

■     i    1     ; 

M        ;        ■ 

""Sp- 

•     ■■     :     . 

1 

: 

'     ■ 

i     ^     i     ; 

'     '     '     ■ 

'•''>, 

1 

\ 

i 

i\ 

V 

:    i    i 

\      , 

1    :    : 

1     ]     ' 

1       ,       ; 

:      ■      1 

i    ,    '    1 

A  i 

r^ 

\ 

1      !      ' 

\ 

f        I 

j\\- 

^\ 



. 





'        ;        I        ■ 

I  '       / 

,       , 

■        '        1 

t-- 

Mr 

Sr-— 

-—- — 

— : 

— —    --■ 



—   J 

\   ■   ■  I 

1    \ 



r  "^  V 



_^_ 

__^-__ 

■-  — 

~~~r 

-— 

/     1     1     \ 

:    1 

rs^ 



I 

/           '     * 

>^ 

1 

/ 

--I 

— __.,. 

— , 

— \ 

1  ._        1  '    _ 



-    -■ 

-p 

—  t  1 

P       r    1 

—- — — 

1     ;     ■ 

i      ■      pv    ■ 

■ 

-  _^ 

1 

h=t 

L_L_.: 

N^ 

\ 

>^,^^ 

/ 

1 /^ 

-1^ 

\ 

-    ~i^s: 

._. 

/ 

/    / 

\ 

_  / 

j     ' 

Li 

V 

.^j — _.. 

r  7     1    r   I 

"       "      ■ 

F 

h^-^^J 

-— : 1 

_ — -^i 



1 

-    Q    =    10.0 

2 

-    Q    =       4.0 

3 

-    Q    =       1.0 

Fig.  5.25  -  BPI  Amplitude  Variation 
(R=1.6k) ,  f=3.83khz 


132 


1    !   i  ■' 

:    :    1    1 

'        '       •        ' 

1  ■  ' 

'   ■ 

■ 

:    ■ 

.    y      .         .,     . 

-   ' 

.     ,     ; 

1       '    1 

;      1      . 

1     1 

1 

1    '    1 

i    !    1    • 

'  i  1 

:      1      !      ! 

III 

. 

2dblai 

i    ; 

] 

ill! 

1  !  ' 

1  i  I 

i    ;    ' 

1    i    ■    : 

•    i    I    ! 

1    ■       i 

1  i  1  i 

■  III 

1    ■    1    ; 

1    '    1 

1 

1         i          ' 

i     i     !     ■ 

:      1      ' 

i      1      ; 

1    !    '    ' 

ill' 

1      I 

1 

i    ' 

!     ■     j 

■       : 

1     '■     \ 

;             ;             , 

1 

1    i       ' 

i  j  1 

!  i    i 

1     1     1 

j    i    1    1 

:     '     1 

1       ■       1 

•     1     '     : 

■     !     1 

1 

1  1      ; 

'     1     ! 

1      i 

i    :    ; 

1    1    ;    i 

i 

!  '  '   1 

■     I     .     ' 

]      1 

1 

1    '    ' 

■  1      ' 

1  1  1 

i  1 

\    \    ' 

1 

; 

■    !    1    ; 

'     :     1 

'     ■     1 

'11. 

III! 

1     !     1     1 

:     ■      : 

!       1       1 

1     ,          : 

'  '  1 

I    ■    !    i 

1     1     i 

'     1     i     1 

1  ! 

!      !      ■ 

;     ■     1 

;  '  '  1 

i    '    i    1 

1     1 

1       !       1 

,     '     ' 

—- — 

1    1    !    ! 

III! 

i     1 

1 

1       ;       i 

1 

r^— • 

^v^-^ 

1 

Mil 

,     ; 

^^ 

.        , 

'  1  .  ! 

^ 

>N,_ 

^.._ 

--—  - 

/ 

^ 

"v^ - 

_ 

! 

X 

■    1 

/  '■ 

J 

r__ 

^"^^ 

/ 

>v 

^*s^ 

/ 

LTirr 

\^- 

^_- 

-.-  ' 

— ■ 

"    *"                'Yi     :     ■ 

•    ■    \    ' 

1       ' 

"t 

■    ■    i 

1 

I 

X 

!    1 

/ 

^y^            -^^ 

^, 

/ 

/       ^ 

\ 

1 

/ 

/'- 

_ 

r    / — 

/    -    -      d 

\. 

/       J 

—  -  Xj \    — 

/ 

/       / 

\ 

1 

/     /  _ 

, .  .     V.  -^.^ 

\         i 

J          1 

t 

-  /  ' 

'     " 

\-     4 
\  ..  1    .... 

i 

' 

/ 

—     / 

~  1 

7 

.  - 

U^-u\^ 

>^ 

-- 

1  - 

2  - 

3  - 

-  Q   =    10.0 

-  Q   =      4.0 

-  Q   =      1.0 

Fig. 

5.26 

BPF  Amplitude  Response 
(R=1.6k) 


133 


1    !    t    ! 

— 1 — • — - — : — ' 

■    ,         : 

.  .  ■  \ 

•-  ;  - 

1    ■ 

.  ,  1 

,  ,  , 

— T    — —    1           —' 

1 

1   ! 

1    , 

1 

1    '    '    1 

ill' 

■    !    i 

1      !      . 

1  ; 

ill: 

1    1    !    1 

,    •    1    ■ 

■      1      i      ' 

^atrial 

1  »     .  ..:.,  1 

1 

1 

i       i              i 

i  i  M 

i      1      i      ' 

1 

1 

!  1    i 

1  i  1  j 

■    1    i    i 

I   1  '  ; 

.  1  i 

j»>/ ; : 

i    I 

•    '    i 

' 

1      1 

1     :     i     ' 

i  ■  '  i 

1   ■ 

1  ;  i  ' 

1    1   ! 

i  i  i  i 

i 

1  ' 

1        ]        ; 

!      i      ' 

i    '    '    ' 

i      '      ; 

!    1    i    i 

;  1 

>     '     I 

i   1   ; 

i  M  : 

Mm 

i  i  i  i 

1 

i 

I        i        1 

1  ;  !  : 

I'll 

'        1        • 

1   ,   , 

'    . 

— 

_^-^--^ 

i   1   ' 

!    1 

3:       i 

:    1    ■ 

;     1     ' 

.              . 

:    1 

__.__ — 

;      1 

^^^ 

/    — ~~ 

/"V 

-r-T 

s^'^ -::!:— 

— 

1     !     1 

'      '      i 

1     ' 

1                   J 

f^ 

'  \ 

ill! 

\ 

'      i 

>w        1 

'       \ 

'      ' 

^>v 

\ 

'S^ 



1              / 

i         . 

k,^^-^— 

------ 

1     ;    ^ 

\~" 



1              ,       : 

I     .A 

1 

J ^^N.- - 

' 

\ 

V 

— r— — 

1                      J 

h\~ 

/ 

/ 

1 
I    ■    ■ 

-     r 

/- 

/                 L 



^*" 

..    - 

!    !    ;    ■ 

, 

/       I 

-  \ 

! 

h-\ 

-- 

— '- 

' 

■      : 

-    /- 

f       / 

._.  /       / 

1 

\ 

^_~;._ ::::.-,  ;:;• 

' 

f 

—  J 

'-/- 

t-^ 

:Vr  [-...-...:- 

/ 

/ 

L. . ',    ,  /      / 

\ 

\ 

Q 
Q 
Q 


10 
4 
1 


Fig.  5.27  - 


BPF  Amplitude  Response 
(R=1.6k),  f=9.khz 


134 


';<  F  R 
XFR 

FCl 
FCl 

4  .  • 
C  E  H  T  E  R 

3dB 

'J 

tr 
■- 

dB.' 

0  0  / 

Diy 

r~"- 

.Z_ 

/ 

/ 

,1 

x 

\. 

1 

' 

\ 

"    ~1 

"1. 

i 

'"■■i 

/ 

'> 

I."J. 

•■■.  Gi  H  2 

ti  k:  R  ^   7  8  0  0  H 


Z  5     K  H  z  /•■ 
B  I'J  :  2  0  0  H  2 


1.  Phase  Response 

2.  Amplitude  Response 


:<  F  R 

XFR 

FCl 
FCl 

■  H    M  f. :  R  •■ 

'  U  ■■      0  o 

4  . 
C:  E  H  T  E  R 

3dB 

d  B  y 

0  o  / 

D  I  ',.' 
D  I  '...' 

f---,- 

.-2. 

_-L 

i 

! 

'■■- 

\ 

S 

— »-— 

.  - 

— — 

■-•--■ 

-'--- 

\ 

•-.  0  H  z 

M  K  R  ■■         6  0  0  0Hz 


2  5  K  H  z 
Bli--      2  0  0  H 


Cb) 


Fig.  5.8  -  BPF  Phase/Amplitude 
Response 


135 


X  F  R     F  C  T  W     M  K  R  :  3  .  6  d  B        2     d  B  '■  D  I  './• 

X  F  R     F  C  T  N  =      (do     CENTER  50-^'  y  D  I  V 


r\ 

V 
1'. 
1 

/ 

'\ 

1  - 

\ 

■ 

rf 

"r 

I 

\ 



■••     Q     H  2 

f1  f(  R  :         2     b  CI  0     H 


2  5     f- ■  H  z 
B  I'J  •      2  0  0     H 


1.  Phase  Response. 

2.  Amplitude  Response 


\<  F  R 
X  F  R 

FCl 
FCl 

'  H     M  f:::  R  ■■ 

■  ^.^ :    r-  ^ 

1     . 

C  E  H  T  E  R 

3dB 

e 

d  B . 

0  O  .  ■■ 

D I '.; 
DI  y 

r--.. 

', 

1 

1 
1 

\, 

1 

1   -.-, 

\ 

2. 

.-_'•-  ,-_•-" 

-■-"-■*-*— 

•"-"-"■" 

'"—-—.'  ' 

"■■""" 

V 

1 

1^ 

■ 

••■.  0  H  z 

M K  R  ^    2  6  0  0  Hz 


2*5  KHz 
BN=  2  00  H 


Fig.  5.29  -  BPI  Phase/Amplitude  Response 


136 


^^FR 

XFR 

FCl 
FCl 

'\\  ■        0  0 

2  .  I 
C  E  f-^  T  E  R 

5dB 

2    dB/ 

5  0  o  ,/ 

on; 

r"~ 

'■—- ■ 

...2 

h 

Y 

V 

f 
— f — 

\ 

_,_ 

----- 

-..  ..■^.., •..-_.- 

\ 

-V 



1 

1 

1 

\. 

1 

-.  0  H  z 

M  f(  R  :    ^  8  0  0  Hz 


2  5  KHz  . 
B  l-J  ■■     2  0  0  H : 


1.  Phase  Response. 

2.  Amplitude  Response 


■A  F  R  F  C  T  H  n  K  R  ■■  2  .  3  d  B   2  d  B  ■ '  D  I  '..' 

:■<;  F  R  F  C  T  hJ  =  0  o  C  E  hU  E  R        5  0  o  ..•  d  i  y 


r-- 

..jZ 

/  \, 

1 

\ 

-^-^_.- 

,■^-.-.- 

---—'■- 

..-^.- 

/ 

\ 

1.. 











/ 

1 

/ 

s   0   H  2 

N  K  R  =    6  0  0  0  Hz 


2  5  K  H  z 
B  W  =  2  0  0  H 


5.30  -  BPF  Phase/Amplitude  Response 


137 


X F  R     F  C  T  H     n  K:  P  :      -            .  8  d  B        2     d  B  .'  D  I  '^ 
XFR     F  C  T  H  ^      0  o     c  E  H  T  ER 5  0  0/01  \> 


■■■■■     0     H  2 

M  f-:  R  ■■        2    e  0  0    Hz 


2  5    f:;:  H  z    .  ■■■ 

B  U  ■■      2  0  0     Hz 


1.  Phase  Response. 

2.  Amplitude  Response. 


K  F  R 

>=:fr 


F  C  T  H  M  f :  R  ^  OFF  S  C  R  L  E  2  d  B  ■■■ 
FCTH   oo  CENTER        5  0  o,.- 


D  I  '..' 
D  I  '...' 


4 

:■'     \ 

^r— - 

'"--1 

I 

! 

1 
1, 

1 

\ 

1 
1 

'l 

— — 

_.J^_ 

h-  - 

..-..- V- 

-"■_/\_l*-*- 

/ 



\^ 





1 

y 

\  0 

f'1  k:  r 


Hz 


6  0  0  0  Hz 


2  5  K  H 
B  N  •  2  0  0 


Hz 


Fig.  5.31  -  BPF  Phase/Amplitude  Response 


138 


;^FR 
XFR 

FCl 
FCl 

'W  ■■      0  o 

OFF  S  ( 
CENTER 

::nLE 

'J 

rj  B .  ■ 

0  O  y 

D  I '.; 
D  I '.; 

/  \^ 

\ 

r 

■-■-f 

'    1 

I 

\ 

1 

\ 

\ 

i 

1 

_  :z.. 

.-\.\., 

.-■v^— 

.,...,(, 

-V-' 

% 

/ 

\ 

/ 

\. 

\    i 

0  Hz 
ri f :.  R  :    6  0  0  0  Hz 


2  5  f- •  H : 
B  N  ■■      2  0  0 


Hz 


1.  Phase  Response. 

2.  Amplitude  Response. 


:•<  F  R 
;•<  F  R 

F  C  T  H  M  K  R  : 

F  r:  T  W  ^      0  o 

OFF  G  ( 
C  E  H  T  E  R 

:  R  L  E  2 

d  B . 

0  O 

D  I  '..' 
D  I  'J 

1 

/ 
1 

\<i. 

1 

f'-" 

-> 

1 
1 

i' 

1 
1 

\ 

-^ 

"i 

1 
\ 

1 

'■■■ 

I 

't 

I 

i. 

\ 

r 

■  "" 

" 

/ 

1 

■•,_ 

/ 

-  -"r-j 

• .  0  H  z 

r'1  K  R  ■         7     6  0  0  H : 


2 5  f(  Hz  / 
B  W  =  2  0  0  Hz 


Fig.  5.32  -  BPF  Phase/Amplitude  Response 


139 


■    1    '    ' 

1    '    i 

t 

,        ;         , 

i   i   1 

ill. 

!      i      ■      1 

'  M  ; 

1   i 

'   i 

I                 .        1 

.500  h 
2  mV  / 

1 

i  ' 

M    !    i 

MM 

:        i 

,     1     ;     , 

!  1  ' 

1      i 

1 

' 

■    1    ; 

kU 

1 

1  ' 

; 

i    '    1    ' 

1  '  ' 

1        !        1 

'    1    i    ; 

I      :      i      ■ 

1    1    ■ 

: 

1     !  ' 

i  M  1 

1        1        <        , 

III 

1        . 

•     1     1     ; 

1  !  h 

;   '   i   • 

1 

1  1  i  1 

i     '     '     i 

1    '    !    ! 

i         !        ■         ■ 

'    1    i 

■   1   i 

i  i  i 

i   ! 

.      '     ,     I 

'     !     ;     ' 

'         I     ' 

' 

1 

i   '   i 

1 

----- — 

i 

1  i  i 

;     1     1 

Mil 

!     i 

'     i 

\ 

;    ;    , 

[ 

i     ; 

'      ;      ■ 

i 

; 

\     ;     ■     • 

1,  _ 

r== 

1     .     !     ' 

^^^ 

• 

\ 

XI 

\ 

N 

I 

> 

I 

3^ 

\ 

2 

Y 

V.    '   ■ 

! 

i\ 

\ 

\               ____.. 

_^  A — 

_i_-^„ 

-\— :. 

-„._.aj^_ 

J  _L'__' _ 

1  r   . 

- 

( — ' 

1    ;    :    . 

— V- 

-  y^^ 

-   V- 

^ 

——r:.:A 

-^-"'L^ 

■   -   ■ 

_[  1    ,    1 

\^ 

X. 

---^x-^ 

^~' 

^ 

\ 

— -T-- — 

— :^--- 

;     ; 

_.. 

— 

'  t  i 

! 

!  '■ 

,     ; 

,'•'.• 

1   ;   1 

'    '    i • 

\ 

1 

_ 

Fig.  5.33  -  Notch  Amplitude  Response 
(R=1.6k) ,  Q=4 


140 


. 

1 

t     1     .     . 

r~            ,       1    -    V.        ^ 

•    1 

!  1  '  ■ 

■    : 

i    J    '    ! 

;     i 

,  1 

\\d\^ 

III! 

1  i  '  1 

i  1 

'    '       1 

j    i    1    1 

(     1 

1  ' 

♦  A^  , 

'        i        '         ' 

1     :     ■ 

•    1       1 

!   !   1   i 

/dt/ 

'   '   ;   ■ 

1     i 

' 

1   i       ' 

:    1    i    • 

;     '     [     , 

'               .        . 

1    1    1 

' 

I'M 

1     1     •■ 

1 

, 

; 

1    j    ' 

■ 

^-— 

1 

'  :  1 

y"v 

, 



_-.. 

^^<^ 

V^'^N  ' 

[\   : 

;      '      ■ 

■ 

""4"" 

^-^    1    \'    '    ' 

'^— - 

iX^z 

: 

' 

, 

1     [^ 

\      \ 

\ 

— 

—  1 

-    - 

( 

I 

\ 



— 

\ 

4- 

r  zV 

^r 

. 

— — 

_._ . 

■ 

1      ■      ^      . 

r 

\ 

—  -  -- 

■— --  • 

,     ■     ; 

1 

__^__       I 

— — ■ — 

^ 

1 

/^      I 

^--j \ 

■    ~ 

'■:::zj^^. 

i_       \ 

^^^^^"^     j 

t-    v 

-     "   4     ^\ 

L 

— - — - 

— -  - 

^-    r 

F— t/^' 

I 

L 

y'^             . 

\ 

/^ 



1  ' 

/    ' 

,  \ 

/i                \ 

_^_\ 

/_            iV 

1  ' 

.... 

\->>1 

K— 

I 

i 

lA-Z, 

---■ 

\i 



'  1 

' 



t 1 

,  1  , 

_        .     ^^ 

.J 

V- 

0.99khz 
1.99khz 
3.83khz 


Fig.  5.34  -  Notch  Amplitude  Response 
(R=1.6k) 


141 


illustrate  the  amplitude  response  for  different  frequencies 
and  for  Q=4,  while  Figs.  (5.35)  and  (5.36)  illustrate  the 
amplitude  response  for  f=lKHZ  and  a  variety  of  Qs.  Fig. 
(5.37)  illustrate  the  phase  response  in  addition  to  the 
amplitude  one. 

5.   All  Pass  (AP)  Realization 

The   topology-control   bitword   100   realizes   the 
All-Pass  filter.   Figs  (5.38),  (5.39)  and  (5.40)  illustrate 
the  plotted  amplitude  and  phase  responses. 
B.   CONCLUSION 

The  constructed  circuit  performed  as  predicted  by  the 
theoretical  analysis  and  the  computer  simulations.  This 
means  that  it  realized  all  the  desired  filtering  transfer 
functions  LP,  HP,  BP,  N,  and  AP.  The  effect  of  interference 
of  the  control  switches  nonideal  performances  which  is  more 
severe  at  high  Qps  can  be  minimized  by  increasing  the  values 
of  R  (meaning  that  the  values  of  Rq  also  increases  and  the 
values  of  Cs  decreases). 
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Fig.  5.38  -  APF  Amplitude  Response 
(R=1.6k) ,  f=2.8khz 
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VI.   COMBINING  HIGHER  ORDER  SECTIONS 

By  cascading  two  or  more  programmable  filters,  higher 
order  transfer  functions  can  be  obtained.  Fig.  (6.1) 
illustrates  two  cascaded  GIC  programmable  filters. 


"J,. 


/ai_L 


Fig.  (6.1)  Two  Cascaded  GIC  Programmable  Filters. 
(Each  Black  box  Stands  for  the  Network 
of  Fig.  (3.10) 


25  different  combinations  of  the  two  individual  transfer 
functions  Ti(s)  and  Tj(s)  can  obtained  as  it  is  indicated  in 
Table  (4.3). 

For  Lp-LP  combination  and  for  ideal  theoretical  case 
(Al=A2-->oo)  a  fourth  order  low  pass  filter  can  be  obtained 
with  transfer  function: 


(6.1) 
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The  computer  simulation  of  the  fourth  order  transfer 
function  of  a  nonideal  theoretical  filter  vs.  the  second 
order  case  is  illustrated  at  Fig.  (6.2).  The  experimental 
results  taken  from  the  properly  designed  and  built  circuit  as 
illustrated  in  Fig.  (6.3)  are  shown  in  Fig.  (6.4). 

Using  the  same  procedure  as  above  a  fourth  order  HP-HP 
combination  for  the  ideal  theoretical  response  is  given  by 


\\i-^  ^  ^ 2 (6.2) 


which  at  Wp  takes  the  complex  value  of 


\iG^V)---^^^f"  (6.3) 


with  magnitude  of 


Fig.  (6.5)  illustrates  the  fourth  order  nonideal 
theoretical  HP  filter  response  vs.  the  second  order  one. 
Both  at  Qp=2  and  f p=8KHZ  ,  while  fig.  (6.6)  illustrates  the 
experimental  responses.   For  Qp=2  according  to  (6.4)  and 
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LRF  AMPL  RESP0NSE(Q=2/P=4K) 


FREQUENCY  HZ 


Fig.  6.2  -  LPF  Fourth  Order  vs.  Second  Order 

Ideal  Response  from  Computer  Simulation 
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H.P.F  AMPLRESP0NSE(Q=2/P=4K) 


FREQUENCY  HZ 


Fig.  6.5  -  "Ideal"  4th  vs.  2nd  Order  HPS  Responses 
From  Computer  Simulation. 
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4th  Order  vs.  2nd  Order  HPF  obtained 
Responses  from  the  Constructed  Circuit 
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(2.4)  a  difference  of  12dB  were  expected  between  the  second 
and  fourth  order  filter  which  is  approximately  the  case  in 
both  computer  simulation  and  experimental  results. 

The  BP-Bp  combination  leads  to  the  theoretical  "ideal" 
fourth  order  transfer  function: 


4  ^  ^^-^   (6.5) 


which  takes  the  following  value  at  ujp: 


HLj  o:)?^r  4  --  A2aV3  (6.6) 

According  to  this,  a  magnitude  response  of  12  dB 
approximately  at  uop,  was  expected  from  both  experimental  and 
computer  simulation  results.  Fig.  (6.7)  illustrates  the 
simulated  response  of  a  fourth  order  nonideal  HP  filter  vs.  a 
second-order  one.  A  difference  of  5.8  dB  instead  of  6  dB  can 
be  observed.  Fig.  (6.8)  illustrates  the  experimental 
response,  where  a  difference  of  5.4  dB  appears  basically  due 
to  the  interference  of  the  control  switches  and  the 
nonideally  matched  values  of  the  capacitors  which  control  the 
fp  selection. 
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B.P.F  AMPL  RESP0NSE(Q=4/P=8K) 


FREQUENCY  HZ 


—  cascade 
single 


Fig.  6.7  -  Fourth  Order  vs.  Second  Order  BPF 
Amplitude  Responses  from  Computer 
Simulation. 
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The  N-N  combination  for  the  ideal  theoretical  case 
(Al=A2-->oo)  results  the  following  transfer  function: 

T-2Ci)=        (^''^tf  (6.7) 

I 

where  cDn  is  the  Notch  frequency. 

Fig.  (6.9)  illustrates  the  fourth  order  "nonideal"  Notch 
filter  simulated  response  vs.  the  second  order  one. 

The  AP-Ap  combination  for  the  theoretical  ideal  case 
result  in  the  following  transfer  function: 


(6.8) 

'? 
which  takes  the  following  values  for  S=0  and  oo: 


TlGb^^i  (6.9) 


with  amplitude  and  phase  of 


iMlib^Ui-.bdV,,  l:[iM_  =  3^o' 


(6.10) 
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N.F  AMPLITUDE  RESP0NSE(Q=2./Fr8K) 
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Fig.  6.9  - 


Fourth  Order  vs.  Second  Order  Notch 
Filter  Amplitude  Response  from  Computer 
Simulation. 


160 


and 


Mjoo^.-al5  (6.11) 


with  amplitude  and  phase  of 

lTl(joQ^l.  7c\b    .    /     ^Uj^co)        ,|fAU  (6.12 


The  computer  simulation  of  nonideal  fourth  order  all  pass 
filter  vs.  a  second  order  one  which  is  illustrating  at  Fig. 
(6.10)  matches  the  above. 

Fig.  (6.11)  illustrates  the  resulting  Chebychev  filter 
from  a  BP-BP  combination  with  Qp=4  and  different  frequencies 
(fpl=8KHZ  fp2=10KHZ),  while  Figs.  (6.12)  and  (6.13) 
illustrate  the  resulting  Chebychev  filters  from  the  designed 
and  built  circuit.  Fig.  (6.14)  illustrates  the  resulting 
response  from  a  HP-LP  combination. 
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A.P  AMPLITUDE  RESP0NSE(Q=2) 


9id(f 


FREQUENCY  HZ 


Fig.  6.10  - 


APF  Fourth  Order  vs.  Second  Order 
Amplitude  Response  from  Computer 
Simulation. 
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B.P.F  AMPUTUDE  RESP0NSE(Q=4) 
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Fig.    6.11    -    Chevycheb    Response    for    f-j^  =  8k. 


f2=10k 
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Fig.  6.12  -  Chevychev  Response  Obtained  by 
Cascading  BP-BP 
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Fig.    6.13    -    BP-BP    Response    Q,=Q^=' 


(f^=3khz,    f2=15khz) 
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'ig.     6.14    -    HP-LP    Response    Q-,=Q2  =  4 


(f^=3khz,    f2=9.0khz) 
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VII.   APPLICATION  OF  THE  PROPOSED  GIC  PROGRAMMABLE  FILTER 
IN  FREQUENCY  HOPPING  SYSTEMS 

A.   BACKGROUND 

1.   General  Description  of  Frequency  Hopping  Signals 

Frequency  hopping  is  a  spread  spectrum  modulation 
technique  used  to  generate  many  possible  carrier  frequencies 
over  a  large  bandwidth.  Of  all  the  possible  carrier 
frequencies,  only  one  is  selected  at  a  given  time.  However, 
all  frequencies  are  eventually  selected  during  some  time 
interval. 

Frequency  hopping  (FH)  may  be  pictured  as  an  RF 
carrier  whose  center  frequency  is  "hopped"  over  many 
frequencies.  The  hopping  may  be  either  in  a  simple  sequence 
or  a  pseudorandom  sequence. 

The  hopping  rate  of  a  frequency  hopping  system  does 
not  affect  the  bandwidth  of  the  output  spectrum.  In  a  direct 
sequence  system  the  chip  rate  determines  the  total  bandwidth. 
In  a  frequency  hopping  system,  however,  the  bandwidth  is 
determined  by  the  highest  and  lowest  frequencies  of  the 
frequency  hopped  carriers.  For  example,  if  the  highest 
frequency  carrier  is  at  15  MHz  and  the  lowest  frequency 
carrier  is  at  10  MHz,  the  total  signal  bandwidth  is  5  MHz. 
This  is  the  bandwidth  regardless  of  the  hopping  rate.  This 
allows  wideband  spread  spectrum  signal  generation  at  low 
hopping  rates. 
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2.   Signal  Generation 

Frequency-hopped  signals  may  be  generated  in  several 
ways.   The  different  methods  are  classified  into  two  groups: 

(1)  Direct  synthesis,  and 

(2)  Indirect  synthesis. 

One  important  aspect  of  frequency  hopping  synthesis 
is  coherency.  coherent  signal  synthesis  is  defined  as  the 
establishment  of  a  known  and  repeatable  phase  each  time  a  new 
frequency  is  hopped  to.  Non-coherent  signal  synthesis  is 
defined  as  the  establishment  of  a  random  or  unknown  phase 
each  time  a  new  frequency  is  hopped  to.  some  techniques, 
direct  or  indirect,  of  signal  generation  can  be  used  as  a 
coherent  frequency  source.  In  other  techniques,  the  changing 
of  frequencies  creates  non-coherent  sources. 

If  a  frequency  hop  system  is  a  coherent,  it  will  have 
a  signal-to-noise  advantage  over  a  non-coherent  system. 

a.   Direct  Synthesis 

The  direct  approach  to  signal  synthesis  utilizes 
techniques  which  enable  direct  synthesis  of  different 
frequencies.   Examples  of  direct  synthesis  techniques  are: 

(1)  Frequency  mixing,  and 

(2)  Surface  acoustic  wave  devices. 

Frequency  mixing  for  single  synthesis  is  a  common 
technique  used  to  generate  many  different  frequencies.  An 
example  of  the  frequency  mixing  technique  is  show  in  Fig. 
(7.2). 
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FlGURr?.!  BLOCK  DIAGRAM  OF  A  FREQUENCY  HOP  MODEM 

A  pseudorandom  code  generator  selects  one  of  many  transmit  frequencies  during  a  sma H 
till!"  interval.  Tlie  traffic  "lodulatos  t'le  carrier  frequency  v/hich  is  spread  over  many 
different  frequencies  hy  the  horp'ng  action.  The  receiver  dehops  the  input  signal 
into  a  narrowband  If.  Tlie  code  synchronizer  locks  the  pseudorandom  code  generator  in 
the  receiver  to  the  received  signal.  A  data  demodulator  removes  the  traffic  from  the 
ir  ami'l  i  f  ier  output . 
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In  Fig.  (7.2),  an  RF  switch  selects  one  of 
several  frequency  inputs.  Two  of  these  input  signals  of 
different  frequencies  are  multiplied  together  to  generate  a 
new  output  frequency.  The  device  used  to  multiply  the  two 
signals  together  is  called  a  frequency  mixer.  when  two 
frequencies  are  mixed,  the  sum  and  the  difference  of  the 
frequencies  are  generated.  In  order  to  select  only  one  of 
these  frequencies,  a  "filter  is  used  to  reject  the  unwanted 
frequency."  A  filter  tuned  to  the  desired  frequency  would 
allow  selection  of  that  frequency  while  rejecting  the  other. 
By  selecting  the  mixing  frequencies  in  the  proper  order,  the 
output  frequency  can  be  stepped  through  several  different 
frequencies.  At  each  frequency  mixer  output,  a  filters  is 
required  to  reject  unwanted  frequencies.  The  filters  may 
require  a  short  time  for  the  signal  to  stabilize  after  it  is 
selected.  The  time  required  for  the  filter  to  stabilize  at 
each  new  frequency  may  ultimately  determine  the  maximum 
hopping  rate  of  the  direct  frequency  synthesizer, 
b.   Indirect  Synthesis 

The  indirect  method  of  signal  synthesis  is 
defined  as  frequency  synthesis  through  the  use  of 
phase-locked  oscillators.  One  common  indirect  method  for 
synthesis  is  shown  in  Fig.  (7.3). 

In  this  circuit  a  phase-locked  loop  is  used  to 
generate  the  numerous  carrier  frequencies.  The  phase-locked 
loop  has  an  internal  oscillator  whose  output  frequency,  Fo, 
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r  I  CURE  7.2  DIAGRAN  OF  A  DIRECT  FREQUENCY  SYNTHESIZER 

A  piciidoramlum  spqupnce   <;clpct";    diffprent   combinations   of   F^    througli    F    ,   which  are  mixed 
toqpther   to    form  a   nnw    frenwpncy.      Each  new    frequency    Is   mixed  with   the   rf  oscillator 
for   tran-,miss1oii  as  one  of   the  matiy   hopping   freguencles . 
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FIGURE  7.3   INDIRECT  SIGNAL  SYNTHESIS  USING  A  PHASE-LOCKED  LOOP 

When  the  reference  frequency,  F , ,  Is  at  the  same  frequency  as  F, *,  the  voltage- 
controlled  oscillator  produces  a  constant  output  frequency,  Fq.   Sjnce  Fq  Is 
divided  by  the  programmable  divider,  Fg  Is  equal  to  N  times  F,  .   By  changing 
the  divider  ratio,  N,  many  output  frequencies  are  possible. 
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is  shown  in  Figure  (7.3)  as  the  phase-locked  loop  output.  The 
divide-by-N  circuit  divides  this  oscillator  frequency  by  a 
selected  number,  N. 

The  phase-locked  loop  internally  adjust  its 
output  frequency  Fo  so  that  Fl*  is  the  same  frequency  as  the 
reference  frequency  Fl .  If  the  divide-by-N  circuit  output 
frequency  is  initially  lower  than  the  reference  frequency  Fl, 
the  oscillator  output  frequency  is  automatically  increased 
until  Fl  and  Fl*  are  identical.  When  this  occurs,  the 
oscillator  output  frequency  will  become  stable  and  remain  at 
that  frequency  until  the  number,  N,  changes.  When  this 
number  changes,  the  oscillator  frequency  is  again 
automatically  adjusted  so  that  Fl  and  Fl*  are  again 
identical. 
B.   PROPOSED  USES  OF  PROGRAMMABLE  FILTER 

A  wide  field  of  applications  exist  in  FH  systems  for  the 
programmable  GIC  filter.  The  outstanding  performance  of  the 
filter  (including  sensitivity  and  stability)  and  its  high 
speed  response  to  the  different  inputs  (due  to  the  use  of 
CMOS  integrated  circuits)  make  it  very  exceptional  in  this 
field.  The  first  proposed  use  is  indicated,  in  Fig.  (7.4). 
The  figure  illustrates  a  receiver  of  a  frequency  hopping 
demodulator.  The  received  frequency  hopped  signal  after 
heterodyned  by  the  RF  mixer  passes  through  the  GIC  filter  in 
a  BP  realization  (at  this  application  the  topology  control 
network  does  not  need  to  exist  since  the  BP  realization  is 
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Fig.  7.4  - 


The  Proposed  Receiver  of  A  Frequency 
Hopping  System  Using  the  Programmable 
GIC  Filter 

174 


the  only  type  to  be  used) .  The  frequency  shift  of  the  filter 
is  controlled  by  the  synchronized  pseudorandom  code 
generator.  This  code  generator  also  controls  Q  through  some 
interfaced  binary  logic  in  order  to  correct  the  amplitude 
reduction  which  appears  at  high  frequencies.  The  frequency 
shift  problems  can  be  easily  corrected  by  the  use  of 
composate  Op.  Amp.  Fig.  (7.5)  as  it  is  extensively  analyzed 
in  Refs.  [32],  [34],  and  [35], 

The  programmable  filter  can  also  be  used  in  the  direct 
frequency  synthesizer  as  it  is  illustrated  at  Fig.  (7.2).  At 
each  frequency  mixer  output  exists  the  need  of  a  filter 
required  to  reject  the  unwanted  frequencies.  The  frequencies 
to  pass  are  not  always  the  same,  but  they  hop.  the  BP 
realization  of  the  filter  is  used  which  center  frequency  can 
be  controlled  accordingly.  The  filter  may  require  a  short 
time  for  the  signal  to  stabilize  after  it  is  selected.  The 
time  required  for  the  filter  to  be  stabilized  at  each  new 
frequency  may  ultimately  determine  the  maximum  hopping  rate 
of  the  direct  frequency  synthesizer. 
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VIII.   CONCLUSION 

The  novel  design  described  here  has  resulted  in  a 
universal  prograinmable  filter  than  can  be  digitally 
controlled  to  realize  almost  any  practical  filter 
specifications.  This  is  done  through  the  use  of  CMOS 
switches  controlled  by  binary  codes  to  program  the  order  of 
the  filter,  the  filter  topology,  the  filter  center  frequency 
and  selectivity.  The  design  procedure  required  developing 
optimum  switching  arrangements  for  the  minimum  redundancy  in 
components  and  least  dependence  of  the  filtering  function  on 
switching  imperfections  such  as  switches  stray  capacitances 
and  non-zero  and  nonlinear  switch-on  resistance.  The 
sensitivities  of  Qp,  Wp  are  found  to  be  low  with  respect  to 
the  passive  and  active  elements  variations.  The  experimental 
result  show  close  agreement  between  theory  and  practice. 
Further,  these  results  indicate  that  these  realizations  are 
insensitive  to  temperature  and  power  supply  variations.  A 
wide  field  of  applications  exists  for  the  programmable  filter 
beside  the  one  discussed  in  Chapter  VII. 

(1)  Word  recognition  and  speech  synthesis; 

(2)  Music  applications; 

(3)  Signal  processing  in  communication; 

(4)  Adaptive  balancing. 
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Further  investigation  is  needed  to  develop  a  programmable 
switched  capacitor  realization  that  can  allow  frequency 
scaling  by  changing  clock  frequency.  Work  can  be  also 
extended  for  developing  a  wide  bandwidth  programmable  filter 
using  the  composite  operational  amplifier  technique  proposed 
by  [39].  Such  implementation  would  lead  to  a  very  useful 
monolithic  device  at  moderate  cost. 

The  research  has  yielded  a  paper  that  was  presented  at 
the  19th  Annual  Asilomar  Conference  on  circuits,  systems  and 
computers,  Monterey,  California?  November,  1985  (Appendix  C). 
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APPENDIX    A 

$  JOB 

DI MENS I  ON     AT1S(  100)   ,PT1SI  10  0)   ,AT2S(  lOOJ, PT2St  100  I 
*, AT3S( lOO},PT35(lOO),FATlilOO),FAT2(100i,FAr3ilOOJ 

OI MENS  I ON     ATI  I  lOO), PTll  100)  |AT2<  1 00 )  , PT2 t  1 00 )  i A 13 (  1 00 }  , 
»    PT3C  1001  ,FPTl(  i00>  ,FPT2l  lOO)  ,FPT3<  100) , 
*    FTFl  I  100)  tFTF2llQ0)  , F TF3 I  1  00  )  t FP H  1 0 0 )  , FP2 (  I  001  ,FP3( 100) 

COMPLEX     TlS|T2S|T33|YlS,Y2S,Y3S|Y4S,YSS,Y6StY7S,YaS,DS,KlS,K:2S, 
*K3S,K4S,KSS,K6S 

COMPt_EX     Tl,T2,T3,Yl  ,  Y2  ,  Y3  ,  Y4,  Y5  ,Y6,  Y7  ,  Ye  f  D»  K  1  ,  K2  t  K3  ,K  *  ,  K5  ,  <  6 

01  MENS I ON     FR(   100)  ,3TIS(  100),HTIS<  100)  ,3T2Sl  100)  ,HT2Sl  100) 
*,3T35(  100)  ,HT3S{  100) ,FBTli  I  00 )  |F8T2 (  I  00 )  ,F8T3i  100) 

DI  MENS  I ON     BTl 1  100)  , HTl {  100)  ,3121  1  CO  J , HT2 (  100)  ,9T3(  100) 
*»HT3f 100) ,FHTlt 100) ,FHT21 100) ,FHT3l 100) 

COMPLEX     CIS,C2S ,C3StGl  S, G2S ,G3S, G4S ,G5S| G6S ,G7S,G8S ,00S,X1S ,X2S, 
*X3S,X4S,X5S,X6S 

COMPLEX    ClfC2,C3,Gl ,G2 t G3 , G4, G5tG6, G7 ,Ga ,00, X 1 , X2 , X3 , X4 , X5, X6,5 

OMEGA     =     0. 

DO     20     K     =     I , 100 

OMEGA     =     OMEGA     ♦     125.^20. 

S      =     CMPLX{      0.0     ,      OMEGA     ) 
C»4:«««4c««^«««<rt:^^«««         DATA      «*«**«*«***:;:**:*****  «$**4:****:**45  *******  ***** 


^l      = 

1600. 

CI   = 

aooE-'i 

RQl  = 

320O. 

R2   = 

1600. 

C2   = 

5oe-9 

RQ2  = 

3200. 

R3   = 

1600. 

C3   = 

I 2.5E-9 

RQ3  = 

3200. 

R4    = 

1600. 

C4   = 

3.12E-9 

RQ4  = 

3200. 

A3  =  A4  =  A3S  = 

A4S  = 

1E5 

Al  =  A2  =  AIS  = 

A2S  = 

1E5 

J  =  4 

I  =  4 

L  =  4 

N  =  4 

c 
c 

C  ******<t«**»***:**seLECT  THE  TYPE  OF  THE  FOURTH  FILTER  *************** 
DO  CASE  L 


CASE 

HPF 

GIS 

= 

1  ./R4 

G2S 

•= 

G  IS 

G3S 

= 

S*C4 

GAS 

= 

GIS 

G6S 

s 

GIS 

G5S 

= 

0.0 

G7S 

'  = 

G3S 

Gas 

= 

l./'RQ4 

CASE 

LPF 

GIS 

= 

l./RA 

G2S 

= 

S*CA 

G3S 

= 

G2S>  l./RQA 

G4S 

= 

GIS 

G63 

= 

0.0 

G5S 

= 

GIS 

G7S 

= 

0.0 

G8S 

=  G1S 

CASE 

NOTCH 

GIS 

= 

l./RA 

G2S 

= 

GIS 

G3S 

= 

S«CA 

GAS 

= 

GIS 

G5S 

= 

GIS 

G6S 

= 

0.0 

GTS 

= 

G3S 
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G7S 

=  G3S 

Gas 

=  l./RQA 

1=  NONE  DO 

C           BPF 

GIS 

=  1./R4 

G2S 

=  GIS 

G3S 

=  S*C4 

G4S 

=  GIS 

G5S 

=  0.0 

G6S 

=  GIS 

G7S 

=  t./RQ4 

G8S 

=  S*C4 

END  CASE 

C 

C«««»«»«9«4i«<:«CE4t« 

SELECT  THE  TYPE  OF  THE  THIRD  FILTER  **:>«******* ****S:*:^i 

DO  CASE  N 

CASE 

C            HPF 

Gl 

=  1./R3 

G2 

=  Gl 

G3 

=  S»C3 

G4 

=  Gl 

.  -,;■  -   G6 

=  Gl 

G5 

=  0.0 

G7 

=  G3 

G8 

=  l./'RQ3 

CASE 

C           LPF 

Gl 

=  1./R3 

G2 

=  S»C3 

G3 

=  G2  ♦  l./Rq3 

G^ 

=  Gl 

G6 

=  0.0 

G5 

=  Gl 

Gr 

=  0.0 

Ga 

=  Gl 

CASE 

C           NOTCH 

Gl  = 

s  1./R3 

G2  : 

:  Gl 

G3   : 

=  S*C3 

G4  r  Gl                                                                      i 

GS  : 

=  Gl 

G6  ■ 

=  0. 

G7  : 

s  G3 

G8   : 

=  i.yRQJ 

IF  NONE  DO                                                                         II 

C           BPF 

Gl 

=  1  ./'R3 

G2 

=  Gl 

G3 

=  S«C3 

GA 

=  Gl 

G5 

=  0.0 

G6 

=  Gl 

G7 

=  1.>'RQ3 

G8 

=  S*C3 

END  CASE 

C  *«**<=<i*«***#**<t*SEL£CT  THE  TYPE  OF  THE  SECOND  FILTER  »*****«***:>*«*«    || 

DO  CASE  I 

, 

CASE 

C            HPF 

YIS 

=  1./R2 

Y2S 

=  YIS 

Y3S 

=  S*C2 

Y*S 

=  YIS 

Y6S 

=  YIS 

Y5S 

s  0.0 

Y7S 

=  Y3S 

Yas 

=  1./RQ2                                                            ] 

CASE 

C           LPF 

YIS 

=  I./R2                                                              1 
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Y2S 

s 

S*C2 

Y3S 

= 

Y2S*     1. 

Y4S 

= 

Y13 

Y6S 

= 

0.0 

Y5S 

s 

YIS 

Y7S 

= 

O.O 

Yas 

=  Y1S 

CASE 

c 

NOTCH 

YIS 

= 

1  .•R2 

Y2S 

s 

YIS 

Y3S 

s 

S«C2 

Y4»S 

= 

YIS 

Y5S 

s 

YIS 

Y6S 

s 

0.0 

Y7S 

s 

Y3S 

YSS 

s 

1./RJ2 

IF     NONE     DO 

c 

3PF 

YIS 

■s 

1  ./R2 

Y2S 

= 

YIS 

Y3S 

= 

S*C2 

Y^S 

= 

YIS 

Y5S 

= 

0.0 

Y6S 

= 

YIS 

Y7S 

= 

1  ./Rq2 

Yas 

=  ; 

S*C2 

END 

CASE 

c 

c ************ 

**** 

SELECT     Tl 

OO    CASE 

J 

CASE 

c 

HPF 

Yl 

= 

l./Rl 

Y2 

s 

Y  I 

Y3 

= 

S*Cl 

Y4 

= 

Yl 

Y6 

= 

Yl 

Y5 

= 

0.0 

Y7 

= 

Y3 

YS 

= 

l./RQl 

CASE 

c 

LPF 

Yl 

= 

l./Rl 

Y2 

= 

S*C1 

Y3 

= 

Y2     ♦     I 

Y4 

= 

Y  1 

Y6 

= 

0.0 

YS 

= 

Yl 

Y7 

= 

0.0 

Y8 

= 

Yl 

CASE 

c 

NOTCH 

Yl 

= 

l./Rl 

Y2 

= 

Yl 

Y3 

= 

S*C1 

Y4 

= 

Yl 

Y5 

= 

Yl 

Y6 

= 

0. 

Y7 

= 

Y3 

Y3 

= 

l./RQl 

IF     NOT« 

DO 

c 

BPF 

Yl 

s 

l./Rl 

Y2 

= 

Yl 

Y3 

= 

S*Cl 

Y4 

= 

Yl 

Y5 

s 

0.0 

Y6 

= 

Yl 

Y7 

= 

l./RQl 

Y8 

= 

S*Cl 

END 

CASE 

/RQ2 


TYPE  OF  THE  FIRST  FILTER  **«t*****^*:ft*«**** 


,/RQl 
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c 

FR{K)     =     aM£GA/6«2e 

S     =     CMPLX< 0.0, OMEGA) 

KIS    =    yiS^tl.*      1-/A2S) 

K2S     =     t  Y2S     ♦     YSS     «•     Y6S) 

K3S    =     (Y4S     ♦     YTS     ♦     Y8S} 

K45    =     Y^S     ♦     Y8S 

<SS     =     YIS     ♦     Y3S 

K.6S     =     YSS     «•     Y6S 
C 

Kl     =     Yl*(l.*      1./A2) 

K2     =     ( Y2     ♦     Y5     ♦     Y6) 

K3    =     I Y4     ♦     Y7     ♦     Y8) 

K4     =     Y4     ♦     V8 

K5     =     Yl      ♦     Y3 

K6     =     YS     ♦     Y6 
C 

XIS    =    G1S»(1.*>     1./A4S} 

X2S    =     tG2S     ♦■     G5S     ♦     G6S) 

X33    =     (G4S    ♦     G^S     ♦     G8SJ 

X43     =    G4S     ♦     G8S 

X5S     =    GIS     ♦     G3S 

X6S     =     G5S     ♦    G6S 
C 

XI     =     GI*(l-*      1./A4) 

X2    =     tG2     ♦    G5     ♦     G6) 

X3     =     1 G4     ♦     G7     ♦     G8) 

X4    =     G4     «■     G8 

X5     =     Gl      ♦     G3  '• 

X6     =     G5     ♦     G6 
C 

c  i 

OS     =<<XS*K:2S*<3S)/A1S     ♦     YIS*Y4S*K6S     ♦■  Y3S*K2S^K3S/ I  A  I  S*  A2S1  ♦ 
«  Y2S*Y3S*Y4S/A2S*'(  Y2S/YIS)C.KIS«Y3S*C  Y7S*Y8S  )     ♦    Y3S*K6S*K3S/A  2S        \ 

Tl  S=(  YIS*Y4S*Y5S     ♦^  Y3S*Y5S«  {  Y7S- YBS  )      ♦Y3S* Y7S$<  Y2S*  Y63  )  ♦  I  Y7  5* 
»  IY1S*Y3S»«K2S)/AISJ/DS 

T2S=(  YIS*Y5S*K4S     -Y1S*Y7S*Y6S     ♦  Y3S*YSS*K4SyA2S    -Y3  S*Y6S«Y7  S/'A2  3  I 

*  ♦Y5S5:Y7S*K5S/AIS  ♦  Y2S*Y33*Y7S  1 /DS  ! 
T3S=  (K1S«Y7S*K2S/AIS     ♦K  I  S*Y4  S*Y5S     ♦  Y3S*  Y7S«:K2S/A  IS*  A2S     ♦Y3S*Y4S  > 

*  «Y5S/A2S     ♦     I riS/Y2S)»KlS*Y3S*Y7S     ♦     K6S* Y3S* Y7S/A2S ) /OS  « 

O     =(  K1*K23:K3)/ Al      ♦     Y1*Y4«K6     ♦  Y  3<«K2  *<  3/ 1  Al « A2  )  ♦     Y2*Y3»Y4/A2    ♦ 

*  <Y2/Y1 )*<! *Y3«< Y7*Ya)  ♦  Y3«K6*K3/A2  I 
Tl=(yi«Y4*YS  ♦  Y3«Y5*(  Y7-Ya  J  ♦  Y3»  V7»  I  Y2*  Yb  )  ♦  (  Y7*  |  Y  I  •■  Y3  )  «<2  ) /A  1  ) /D  j 
T2=1Y1»Y5*K4     -Yl*y7»Y6     ♦Y3*Y5*X4/A2    -Y3*y6*Y7/A2     ♦  Y5*Y7«K.5/A2      ♦•  j 

*  Y2*Y3*Y7  ) /O  j 
T3=(  Kl*Y7*kC2>'A  1     ♦«  l«Y4*YS     ♦Y  3^Y7*K2/A1  *  A2     ♦  Y3*Y4*Y5/'A2     ♦1Y2/YI)« 

*  Kl*Y3«Y7     ♦K6*Y3«Y7/A2 J/O 

C  ! 

C  j 

DOS     =f X1S*X2S*X3S) /A3S     ♦     GIS*G4S*X6S     ♦G3S*X2S*X3S/t A3S* A4S } ♦ 

*  G2S*G3S»G43/A4S*(G2S/G13) «XlS«G3S»(G73*Ga3j      ♦     G3S*X6S*X3S/A4S 
CI  S=(GIS<:G4S«GSS     ♦G3S*GSS*  <G7S-G8S)      ♦G3  S*G7S4t  C  G2S*G6S7  ♦!  G7S« 

»  |GIS*G3S)*X2S)/A3SJ/OOS 

C2S=(GlS*G5SOX4S     -G1S«G7S*G6S     ♦G3S*GSS*X4 S/A4S    -G3S*G6S»G7 S/A4 S  I 

*  ♦G53*G7S*X5S/A3S     ♦     G2S*G3Si:G7S     J/OOS 
C3S=tXIS*G7S»X2S/A3S     ♦X1S*G4S*G5S     ♦G3S*G7SX:X2S/A3S«:A4S     ♦G3S*G4S 

*  *G5S/A4S     ♦     (  G1S/G2S)«X1S4«G3S4G7S     ♦     X  6S*G3S*G73/' A4S  ) /DOS 
C 

0D=< X1*X2*X3 J/A3     ♦     G19G4$X6     ♦Ga* X2*X3/ < A3*A4 | ♦     G2*G3eG4/A4    * 

*  {G2/G1 J*Xl *G3«IG7*G8)     ♦     G3*X6*X3/A4 
Cl  =  f  Gl  *G4*G5*GJ*G5*<  G7-Ga  )  ♦G3*G7*  (  G24>G6)>(  G7*|  G  1  ♦G3  )  *X2  J /A3  J /DO 
C2=i  GI  *G5*X4-G1*G7*G6     ♦G3»G5*X4/A4    -G3*G6*G7/A4     ♦G5*G7*X5/A4     ♦ 

*  G2«G3*G7     l/OO 
C3=f XI *G7*X2/A3    ♦Xl*G4*GS     ♦G3«G7*X2/ A3*A4     ♦G3*G4*GS/A4     ♦1G2/G1»« 

*  X1*G3*G7     ♦X6*G3*G7/A4 I/DO 
C 
C 

ATISIK)  =  2O.*ALOGI01CABS(  TIS)  ) 
AT2S(K:)  =  20  •=»ALOG10(CABS(  T23M 
AT3S(K)     s     20  •*AI.0C10|CABS{  T3Sn 
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APTIS 
ART2S 
APT3S 

AI  TIS 
AIT2S 

AIT3S 


REAL (TIS) 
REAL ( T2S ) 
flEAL<T3S ) 

AlMAGt  Tl  S) 
AlMAGt  T2S) 
AIMAG( T3S} 


PT  1  S  i  <)  = 
PT2S<K)  = 
PT3SfK}     = 


ArAN2t  AI Tl S 
ATAN21 AIT2S 
ATAN21 AIT3S 


ARTIS)*  57.325 
ART2S)*  57.325 
ART3S)*    57.325 


ATlfK)  -     20.»ALaG10ICA8S<Tl  )  ) 

AT21IC)  =     20.*ALOG10(CAaStT2)) 

AT3<K)  s     20.»ALOG10<CAeS(T3) ) 

APT!     =  REALI  Tl  ) 

ART2     =  REALl  T2  ) 

A=?T3    =  REAL!  T3  | 

AITI     =  AIMAG{Tl) 

AIT2     =  AIMAGCT2) 

A1T3    =  AiMAG(T3) 


PTHK)  = 
PT2JK)  = 
PT3TKJ     = 

FATHO  = 
FAT2IK»  = 
FAT3(K)     = 


ATAN21 AITI 
ATAN21 A IT2 
ATAN2  t AITJ 

AT14K) 
AT2tK) 
AT3IK) 


,  ARTl )*  57.325 
,  ART2)*  57.325 
t     ART3)«    57.325 

ATI Sl<» 
AT2SIK) 
AT3S(K» 


FPTI (K) 
FDT21K  ) 
FPT3fK j 


=  PT  It  K) 
=  PT2(<I 
=     PT3IK) 


♦  PTlStKl 

♦  PT2St»C| 

♦  PT3S(K) 


3TISIK)  =     2O.«ALOG10{CA8S(  CIS)  ) 

aT25(X)  =     20  .*ALOGlO(CAaS< C25)  > 

3T3S(K)  =     20.*ALQG10(CA8S{C3S) ) 

ARTIS     =  REALCCXSJ 

ART2S    =  REAL(C2SJ 

A^T3S    =  REAHC3S) 

AI  Tl  S    =  AI  MAGt  CI  S) 

AIT2S     =  AIMAGIC2S) 

AIT3S    =  AIMAG(C3S) 


HTIS<<)  =  ATAN2tAITlS 
HT2SIO  =  ATAN2(AIT2S 
MT3S4K)     =     ATAN2<AIT3S 


ARTIS)*  57.325 
ART2S)*  57.325 
ART3S)*     57.325 


3T1(K)  =     20-*ALOG10<CA8S|C 1 ) > 

aT2CK)  =    20.*ALOGlO<CA8StC2) ) 

aT3(K)  =     20.*AL0G10(CA8SlC3)i 

ARTl     =  REALt  CI  ) 

ART2     =  REAL<C2) 

ART3    =  REALIC3) 

AITI     =  AIMAGICI) 

AIT2     =  AIMAG(C2) 

AIT3     =  AIMAG(C3) 


MTKK)  = 
HT2(K)  = 
HT3CK)     = 

FBTl IK) 
F3T2<>C> 
FBT3(K) 


ATAN2I AITI 
ATAN2 ( AIT2 
ATAN2 (AIT3 


(  ARTl )*  57. 325 
t  ART2>*  57.325 
I  ART3J*  57.325 


=  3T1IK) 

=  aT2|K) 
=  BT3<K) 


♦  BTlStK) 

♦  BT2StK) 

♦  8T3SIK) 
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FHTICK)     s     HTllK)      ♦     HTlSlK) 

FHTZIK)     =     HT2<<)      ♦     HT2S(K) 

FHT3(<)     =     HT3IK)      ♦     HT3S(K» 
C 

FTFltK)      -     ATllK)      ♦        ailfKl      ♦     BTlSfKj*     ATISIKJ 

FTF2KI     =     FAT2<KJ     ♦    F3T2iK) 

FTF3(<)     =     FAT31KI     ♦    FBT3<K) 
C 

FPltK)      =     FHTIIK)      ♦     FPTlCK) 

FP2fK)      =     FHT2<K)      ♦     FPT2{Kj 

FP3<K)      =     FMT3IKJ      ♦     FPT3(K) 
C 

WRITE( 6,66)     FRC<  }  ,F AT  I ( K J ,F3T1 IK ) 
*  ,FTFH»C)  ,HTHK),FHTl  (K) 

C  *  FHT3tK) 

20  CONTINUE 

6€>  F0RMAT(6l  IX,F9.  3)  > 

STOP 
END 
*ENTRY 
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APPENDIX    B  .-  ^  : 

ijoe 

DIMENSION     ATISI  100)  |PTISI  100)  ,AT2S(  100l,PT2S<  100) 
,*,AT3S(  I00),PT3S(  lOO],FArH  lO0),FAT2{  100)  ,FAT3ll00). 

01  MENS  ION    ATI  1  100) ,PT1(  100)  ,AT2(  1 00 ) , PT2 I  100),AT3(100), 
»    PT3(  lOO)  ,FPTi(  100)  ,FPT2<  100)  ,FPT3I  100)  , 
*    FTFl  I  100)  ,FTF2( 100  )  ,FTFJ|  100 ) *FP1(  100)  , FP2 ( 1 00 ) , FP3 !  1 00 ) 

COMPLEX     TlS,T2S,T3S,YlS,Y2S,y3S,Y4S , YSS , V6S, Y7S , Y3S  »OSiKlS»K2St 
*<3S,<4S,>C5S,K6S 

COMPLEX     ri,T2,T3,Yl  ,Y2»  Y3,  Y4,  YSf  Y6f  Y7  ,  Y8  ,  D,  Kl  ,  K2  t  <3  t  K4,  K5  ,  IC6 

OI ME MS I  ON    FR(  lOO)  ,3TlStl00),HTlSl  IOC)  ,aT2S{100)iHT2Sll001 
*»BT3S (  100)  ,HT3S|  I  00 ) ,F8TH  100 )  tFBT2|  100)  tFBT3 J  100 ) 

01  MENS  I ON    OTl (1 00) «HTlt  100)  »BT2(  100) , HT2 (  100}  iBT34  100) 
*,HT3f 100) jFHTK 100) ,FHT2( 100) ,FMT3t 100) 

COMPLEX    Cl S,C2S,C3S,GlS,G2S,G3S,G4S,G5S,G65,G7S,G3S,DOS, XI S| X2S , 

^•XT^      w^c      X^^      IfA^ 

COMPLEX    CltC2,C3tGl  ,G2 , G3 ,G4, G5vG6t G7 ,Ga yOOt XI ,X2 »X3 » X4 , X5t X6 vS 
OMEGA    =     O. 
OO     20    K     =     1 , lOO 
OMEGA     =     OMEGA     ♦     125«*20« 
S     =    CMPLXi      0.0     »     OMEGA     ) 
C««:»«««4:«««^:«<c:»««^<[«         DATA     ***^«tc*4E  *<r^*«**«***<t#*«rSi*«!*4t**  ****«***««t«*^ 


Wl     = 


31     = 

1600. 

Cl        = 

2  0OE-9 

RQl     = 

3200. 

R2        = 

1  6O0. 

C2        = 

50E-9 

RQ2    = 

3200. 

R3       = 

1600. 

C3       = 

12-5E-9 

RQ3    = 

3200. 

R4        = 

1600. 

C4        = 

3-12E-9 

RQ4     = 

3200- 

Wl  S     =     ' 

•SET    VALUE 

M 

A3    =     > 

(\4     =     A3S    = 

A4S     - 

=   wi/s 

Ai     X     i 

1^2     =     AIS    = 

A2S     : 

:     WIS/S 

J     =     4 

1     =     4 

L     =     4 

N    =     4 

c 
c 

C  *«»4:* *********** SELECT  THE  TYPE  OF  THE  FOURTH  FILTER  «*«^4:***^*««:*^* 
OO  CASE  L 


CASE 

HPF 

GIS 

= 

1  ./R4 

G2S 

= 

GIS 

G3S 

= 

S*C4 

G4S 

= 

GIS 

G6S 

= 

GIS 

G5S 

= 

0.0 

G7S 

= 

G3S 

Gas 

= 

1  ./RQ4 

CASE 

LPF 

GlS 

= 

1  •/R4 

G2S 

s 

S*C4 

G3S 

= 

G2S*     1 

G4S 

= 

GIS 

G6S 

s 

0.0 

G5S 

= 

GIS 

G7S 

s 

0.0 

Gas 

=  GIS 

CASE 

NOTCH 

GIS 

= 

1./R4 

G2S 

= 

GIS 

G3S 

= 

S*C4 

G4S 

X 

GIS 

G5S 

s 

GIS 

G6S 

= 

0.0 

/RQ4 
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CASE 


GTS  =  G3S 

G8S  =  1./R34 
Ir     NONE  DO 
C           BPF 

GIS  =  1./R4 

G2S  =  GIS 

G3S  =  S*C4 

G4S  =  GIS 

G5S  =  0.0 

G6S  =  GIS 

GTS  =  1./RQ4 

GSS  =SCC4 
END  CASE 
C 

Ctr««««tc«9«4c4:4:«$«^  SELECT  THE  TYPE  OF  THE  THIRD  FILTER  ***«*****Se*^****Se*^ 
00  CASE  N 
CASE 
C            HPF 

Gl  =  1./R3 

G2  =  Gl 

G3  =  S4C3 

G4  =  Gl 

G6  =  Gl 

G5  =  0.0 

G7  =  G3 

G8  =  1./RQ3 

Gl  =  1./R3 

G2  =  S«C3 

G3  =  G2  ♦  1./RCJ3 

G4  =  Gl 

G6  s  0.0 

G5  =  Gl 

G7  =  0.0 

Ga  -  Gl 
CASE 
C           NOTCH 

Gl  =  I.XR3 

G2  =  Gl 

G3  =  S*C3 

G4  s  Gl 

GS  =  Gl 

G6  =  O. 

G7  =  G3 

G8  =   l./RQJ 

IF  NONE  DO 
C           BPF 

Gl  =  I ./H3 

G2  =  Gl 

G3  =  S«C3 

G4  =  Gl 

G5  =  O.O 

G6  =  Gl 

G7  =  1./RQ3 

G8  =S*C3 
ENO  CASE 
C  **«***********<r*SEL£CT  THE  TYPE  OF  THE  SECOND  FILTER  *************** 
DO  CASE  I 


CASE 
HPF 


CASE 
LPF 


YIS  =  I .•R2 

Y2S  =  YIS 

YJS  =  S*C2 

Y4S  =  YIS 

Y6S  =  YIS 

Y5S  =  0.0 

Y7S  =  Y3S 

YaS  =  1./RQ2 


YIS  =  1./R2 
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Y2S  =  S*C2 

Y3S  =  V2S*  l,/RQ2 

Y^S  =  YIS 

Y6S  =  0.0 

Y5S  ~     YIS 

Y7S  =  O.O 

YaS  =Y1S 
CASE 
C           NOTCH 

YIS  =  l./'R2 

Y2S  =  YIS 

Y3S  =  S*C2 

Y4»S  =  YIS 

Y5S  =  YIS 

Y6S  =  0.0 

Y7S  =  Y3S 

YSS  =  1./RJ2 
IF  NONE  DO 
C            3PF 

YIS  =  1 ./R2 

Y2S  =  YIS 

Y3S  =  S«:2 

Y<»S  =  YIS 

Y5S  =0.0 

Y6S  =  YIS 

Y7S  =  l./'RQ2 

Yas  =S*C2 
ENO  CASE 
C 

C**#******«******  SELECT  THE  TYPE  OF  THE  FIRST  FILTER  ******:5*«tt::*^4:****i 
DO  CASE  J 
CASE 
C            HPF 

Yl  =  l./Rl 

Y2  =  Yl 

Y3  =  S<tCl 

Y4  =  YX 

Y6  =  Yl 

Y5  =  0.0 

Y7  =  Y3 

YS  =:  l./RQl 
CASE 
C           LPF 

Yl  =  l./Rl 

Y2  =  S<«Cl 

Y3  =  Y2  ♦  l./RQl 

Y4  =  Y  I 

Y6  =  0.0 

Y5  =  Yl 

Y7  =  0.0 

YS  =  Y  I 
CASE 
C           NOTCH 

Yl  =  l./Rl 

Y2  =  Yl 

Y3  =  S»Cl 

r^    ~    Yl 

Y5  =  Yl 

Y6  =  O- 

Y7  =  Y3 

Y3  =  l./RQl 
IF  NONE  DO 
C  aPF 

Yl  =  l./Rl 

Y2  =  Yl 

Y3  =  SX«C1 

Y4^  =  Y  I 

Y5  =  0.0 

Y6  =  Yl 

Y7  =  l./RQl 

Y8  =S*C1 
ENO  CASE 
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c 


c 
c 


c 

FR<K)     =     0M£GA/6-2e 

S     =     CMPLX* 0.0, OMEGA) 

<1S     =    YIS*!!.*-      1./A2S) 

K2S    =     ( y2S     ♦     YS3    ♦     Y6S) 

K3S     =     (Y4S     ♦     r7S     ♦     YSS) 

K4S     =     Y4S    ♦     Y8S 

K5S     =     YIS     ♦     Y3S 

K6S     =    YSS     ••     Y6S 
C 

Kl     =     Yl*41.>      t«/A2) 

IC2     =     IY2     ♦     Y5    ♦     Y6) 

K3    =     I Y4     ♦     Y7     ♦     Y8J 

K*     =     Y4     ♦     Y8 

K5     =     Yl     ♦     Y3 

K6     =     Y5     ♦     Y8 
C 

XIS    =    GlS^il.*-     1./A4S) 

X2S    =     <G2S     ♦     G5S     ♦     G6S> 

X33    -     (G4S     ♦     G7S     ♦     G8S} 

XA3     =    GAS     ♦     G8S 

XSS     =    GIS     ♦     G3S 

X6S     =     G5S     ♦    G6S 

XI  =  Gl<r(l.*      t./AA) 

X2  =  !G2     ♦    GS     ♦     G6) 

X3  =  (G4     ♦     G7     ♦     G8) 

X4  =  G4     «-     G8 

X5  =  Gl     ♦     G3 

X6  =  G5     ♦     G6 

OS     =4tClS*K2S*<3S)/AlS     ♦     YIS*Y4S*K6S     ♦  Y3S*K2SC'K3Sy' C  A  I  S*  A2S  1  ♦ 
«  Y2S*Y3S*Y4S/A2S*(  Y2S/Y1S><»K1S*Y3S*I  Y7S*YaS  I     ♦    Y3S*K6S*»C3S/'A2 

Tl  S=<  YlSStY4S«YSS     ♦YSS^YSS*  (Y7S-Ya3)      ♦  Y3S«  Y7S*f  Y2S*  Y6S  )  ♦  I  Y7  5* 

*  I  Y1S*Y3S)«K2S)/AISJ/DS 
T2S=(  YlSCY5S<:K4S    -YIS*Y7S*Y6S     ♦  Y3S«Y5S*K4S/ A2S    -Y3  S*Y6S*Y7  SyA23 

»  ♦Y5SCY7S*K5S/A1S     ♦     Y2S*Y33*Y7S     J /DS 

T3S=  IKIS«Y7SX:K2S/AIS     ♦KIS»Y4S*Y5S     ♦  Y3S*  Y7S<:K2S/ A  IS*  A2S     ♦Y3S*Y4S 

*  »Y5S/'A2S     ♦     I  Y13/Y23)«»C1S*Y3S*Y7S     ♦     K6S* Y3S* Y7S/A2S  ) /OS 
C 

O     ={  K1*<2*K3I/Al      ♦     Yl«Y4*K6     ♦  Y  3««K2  5!K  S/f  Al*  A2  )  ♦     Y2*Y3«Y4/A2     ♦ 
«  (Y2/Y1 )*<1«Y3«( Y7*Ya)      ♦    Y3*K6*K3/A2 

T1=(Y1*Y4*Y5     ♦Y3«Yb*l Y7-Ya>     ♦Y3* Y7» ^ Y2>Y6 > ♦ f Y7* | Y 1 ♦ Y3 ) «K2 ) /A  1 ) / 
T2=(Y1#Y5*K4     -Yl<iY7»Y6     ♦  Y3CY  5«K4/A  2    -Y3*Y6*Y7/A2     ♦Y5*Y7«K5/A2     ♦ 

*  Y2*Y3*Y7     ) yO 

T3=i  Kl*Y7<!»C2yA  I     ♦<  1«Y4*Y5     ♦Y  3^Y7*K2/A1*  A2     ♦Y3*Y4*Y5/A2     ♦tY2/Yl} 

*  Kl*Y3«Y7     ♦K&*Y3*Y7/A2 )/0 
C 

C 

OOS     =(  X1S«X2S*X3S) /A3S     ♦     GIS*G4S*X6S     ♦G3S<:X2S*X3Sy' <  A3S*  A4S  1  ♦ 

*  G2S*G3S«G4S/A4S*<G2S/Gl3l *XlS«G3S»(G7S*GaS)      ♦     G3S*X 6S*X3S/A4 
CI  S=  (  GIS<:G4S*GSS     ♦  G3S*G5S*  I  G  7S-G8S  )      ♦G3  S*G7S*  I  G2S*G6S7  ♦  I  G7S* 

*  IG1S«-G3S)*X2S)/A3S)/00S 
C2S=(GIS*G5S*X4S    -G1S»G7S«G6S     ♦G3S*GSS*X4 S/ A4S    -G3S*G6S*G7 S/AASf 

»  ♦G53*G7S*X5S/A3S     ♦     G2S<!G3S*G7S     )/OOS 

C3S=( XXS*G7S»X2S>'A3S     ♦X1S»G4S*G5S     ♦G3S»G7S*X2S/A3S*A4S     ♦G3S*G4S 

*  *G5S/A4S     ♦     f  GlS/G2S)*XlS4tG3S*G7S     ♦     X  6S*G3S*G73/' A4S  >/DOS 
C 

00  =  t  XI*X2*X3|/A3     ♦     G13:G4$X6     ♦G3*  X2*X3/ (  A3*A4  |  ♦     G2*G3«G4/A4    ♦ 
»  (G2/G1 |«Xl «G3*JG7*Ga)      ♦     G3*X6*X3/A4 

Cl=( Gl *G4«G5*GJ*G5*tG7-Ga)^G3«G7*«G2*G6)*( G7«| G 1 ♦G J ) ^X2 I /A3 J /OD 
C2=<  Gl  «G5*X4-G1*G7»G6     ♦G3«(G5«X4/' A4    -G3*G6«G7/A4     ♦G5*G7*X5/A4     ♦ 

*  G2«G3«G7     I /OO 
C3  =  <  X1*G7<!X2/A3    ♦Xl*G4«G5     ♦G3«G7*X2/ A3*A4     ♦G3«G4«GS/A4     ♦(G2/G1|I 

«  X1*63*G7     ♦X6*G3*G7/A4 I/DO 

C 
C 

ATlSfO     =     2O.«Al.OG10f  CA8S(  TIS)  ) 

AT2S(tC|     =     20  •^ALOGlOtCABSi  T2S)  } 

AT3S(<)     =     20.«ALOC10(CAas(  T3S}  } 
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APTIS  =  REAL(T1S| 
ART23  =  R£AL(T2S) 
ART3S    =     ftEAL(T3S) 

AI  TIS  =  A(MAG(  Tl  S) 
AIT2S  =  AIMAG<T2S) 
AIT3S    -     AIMAGC  T3S) 

PTlSiO  =  ArAN2<AITlS  ,  APTIS)*  57.325 
PT2S<K)  s  ATAN2JAIT2S  »  ART2S)*  57.325 
PT3SIK)     =     ATAN2<AIT3S     «     ART3S)«    57.325 

ATKK]  =  20.«AI_OG101CA8S(T1  }  ) 
AT2flCJ  =  2O.*ALOG10(  CAaS(T2)  J 
AT3tK)     =     20.»Al.0G»0<CAaS(T3  J  ) 

ARTI  s  REACI Tl ) 
ART2  =  REAl_(T2) 
ART3     =    REACC  T3 ) 

AITl  =  AI>4AG(TI) 
AIT2  =  AIMAGJT2) 
AXT3     =     AiMAG(T3) 

PTlfK)  =  ATAN2IAIT1  ,  ARTt)»  57.325 
PT2IK)  =  ATAN21AIT2  ,  ART2)*  57.325 
PT3TIC1     =    ATAN2tAITJ     >     ART3)«    57.325 

FATXtK)  =  ATlfKl  ♦  ATIS(K) 
FAT2!K»  =  AT2tKj  ♦  AT2SIK) 
FAT3(K|     s     AT3<K)      ♦     AT3StK| 

FPT1<K)  =  PTHK)  ♦  PTlStKl 
FPT21KI  =  PT2fKl  ♦  PT2St)C) 
FPT3(<)     =     PT3IK)     ♦     PT3S4K) 

ariSIK)  =  20.:»ALOG10(CA8S(  CIS)  1 
aT2S(K)  =  20 .*ALOGlO(CAaS( C2S) ) 
aT3SCK)     =     20.*ALOGIO{CA8S{C3S} ) 

ARTIS  =  REAL(CIS) 
ART2S  =  REAL(C2S) 
ART3S  =  REAL(C3S) 

AI  Tl  S  =  AI  MAG!  CI  SJ 
AIT2S  =  AIMAGtC2S) 
AIT3S    =     AIMAG(C3S) 

HTISK)  =  ATAN2fAITlS  »  ARTIS)*  57.325 
HT2S(<)  =  ATAN2(AIT2S  «  ART23)*  57.325 
HT3S4K)     =     ATAN2IAIT3S     f     ART3S)*     57.325 

3T1(K)  =  20.*At_OGtOf  CABSIC  1  )} 
BT2tK)  =  20.*AL0Gl0lCABStC2) ) 
aT3(K)      =     20.*AL0G10(CABStC3n 

ARTl  =  real! CI  I 
ART2  =  REALfC2) 
ART3    =    REALIC3) 

AITl  =  AIMAGiCll 
AIT2  =  AIMAG(C2) 
AIT3    =     AIMAG(C3) 

HTKKI  ~  ATAN2<AITl  ,  ART1)«  57.325 
MT2(KJ  =  ATAN2(AIT2  »  ART2)*  57.325 
HT3«K)     =     ATAN2!AIT3     ,     ART3J*    57.325 

FBTICKI  =  3T1IK.I  ♦  BTlSfK) 
FaT2<K>  =  aT2(K)  ♦  BT2StK) 
FeT3{<)     =     eT3<K|      ♦     BT3S(K} 
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FHTIIK)  =  HTl(K)   ♦  HTISCK) 

FMT21K)     =     HT2<K)      ♦     HT2SIK) 

FHT3iK)     =     HT31K)      ♦     HT3SlK) 
C 

FTFltK)      =     ATl(K)      ♦        STICK)      ♦     BTlSfK)*     ATlS(Kl 

FTF2C<I     =     FAT2KJ     ♦    F3T2(K) 

FTF3(<J     =     FAT3iKJ     ♦    FBT3CK) 
C 

FPltK)      =     FHTllK)      ♦     FPTl(K) 

FP2<K»      =     FHT24K)      ♦     FPT2{Kj 

FP3<K)      =     FHT31K)      ♦     FPT3tKj 
C 

WRITE  (6,66)     FR(K),FAT1CICJ,F3T1<K) 
*  jFTFllK) ,HT1|K),FHTI (K| 

C  *  FHT3IK) 

20  CONTINUE 

66  FORMAT! 6( 1X,F9« 3) > 

STOP 
END 
GENTRY 
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ABSTRACT 

In  this  contribution  a  general  purpoao 
digitally  controlled  analog  filter  la  presented. 
The  novel  design  I3  a  cascade  of  second  order 
sections  that  are  Individually  programmed  to 
achieve  any  filtering  topologies.  Two-Olnary 
words  are  used  to  control  the  pole  frequency  up 
and  selectivity  3p  of  each  section  Independently. 
Each  second  order  section  Is  a  Generallzed- 
Immlttance  converter  (QIC)  blquads  which  are 
known  for  their  high  stability  and  low  active 
and  passive  sensitivity.  CMOS  switches  are  used 
to  electronically  relocate  the  minimum  number  of 
passive  elements  to  achieve  function 
pr ogr ammab  1  1  1  t y .  Switches  are  also  used  to 
select  the  number  of  cascaded  sections  to 
realize  higher  order   transfer    functions. 
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INTRODUCTION 


The  availability  of  an  analog  filter  with 
digitally  controlled  "programmable"  coefficients 
has  been  the  goal  of  many  researchers  due  to  Its 
several  attractions.  One  possibility  of  a  com- 
pact, versatile  analog  filter  under  remote 
control  opens  up  many  novel  and  Independent 
application  areas.  Also,  when  a  programmable 
fllte-  Is  combined  with  a  permanent  reference 
memory  which  is  user-programmable,  this  would 
form  an  econonlcal  and  versatile  device  for 
dedicated  stand-alone  applications.  The  need 
for  such  a  device  was  motivated  by  advancement 
In  film  and  se-nl  conductor  technologies  as  well 
as  the  continuous  upgrading  of  systems  specifica- 
tions to  t(K«  advantaga  of  tna  avallabls 
technologies   to   the   limits. 

Linear  analog  filtering  finds  many  appllca^ 
tlons,  such  as  speech  processing  (recognition  or 
syntheslsj,  geology.  Instrumentation,  communica- 
tions, process  control,  adaptive  balancing,  etc. 
There  has  been  much  emphasis  on  performing  the 
filtering  function  digitally,  largely  because  of 
the  ease  of  varying  and  optimizing  the  transfer 
function.  However,  for  many  reasons,  such  as 
cost,  size,  signal  processing  complexity,  and 
bandwidth.  It  would  be  desirable  to  perform  the 
filter  function  with  linear  components  yet 
retain  the  flexablllty  of  varying  the  filter 
parameters  digitally. 

Recently,  the  advantages  of  combining 
linear  components    (operational   amplifiers    (OAS), 


resistors  and  capacitors)  and  nonlinear  elements 
(switches)  have  been  demonstrated  using  switched 
capacitor  techniques  [i-3].  In  this  contribu- 
tion, we  are  presenting  the  results  of  realizing 
a  continuous  active  device  using  linear  elements 
and  switches  controlled  by  digital  signals  to 
achieve  a  fully  programmable  filter  [1]. 
Several  programming  features  of  the  proposed 
filter  are  reported.  Tne  first  feature  is  the 
ability  of  the  network  to  realize  the  most 
common  filtering  functions  (function  prograimr.a- 
blllty)  namely:  Low  Pass  (LP),  Band  Pass  (BP), 
High  Pass  (HP),  All  Pass  (AP)  and  Notch(N) 
functions,  using  the  minimal  set  of  elements. 
Tne  second  feature  is  the  ability  of  the  network 
to  progra.li  (Independently)  the  key  parameters  of 
the  filtering  function  chosen  (parameter  program- 
mablllty;  namely;  the  pole  resonent  frequency 
(upi  and  selectivity  (ipl.  Finally  the  ability 
to  program  the  network  to  cascade  several  sec- 
tions to  achieve  higher  order  filter.  All  of 
the  above  programmabl 1 1 ty  features  are  performed 
Independently  to  realize  a  universal  filtering 
network . 

2.       DtSICN    ANAl.YCIS   OF   THE    PROPOSED   FILTERS 

The  basic  active  network  considered. here  as 
the  heart  of  the  programmable  filter  Is  the 
second  order  Generalized  Immlttance  Converter 
(OIC)  structure  TS],  Fig.  1,  wnose  superior 
performance  was  establlsned  In  tne  literature 
[6].  Tlie  general  transfer  function  realized 
T( s)    is   gl ven   by: 

T(3)  •■  H(s  )/Dl3  >-(ao->a,3»a32)/(b(3»b,3«b2s2)  (1) 
The  CIC  transfer  functions  of  Fig.  1  assuming 
non-ideal   OAs  are  given  by 


stability   and  sensitivity   analysis: 

An  Important  criterion  of  a  realization  Is  its 
sensitivity  to  element  variations.  The  QIC 
sensitivity  analysis  has  shown  to  be  as  good  or 
better  than  all  competitive  second  order  net- 
worKs  [6].  While  the  QIC  stability  can  easily 
be  demonstrated  since  In  all  the  transfer 
functions  (2),  the  coefficients  of  D(3)  are  seen 
to  remain   positive   for  any  OA  mismatch.      This    Is 
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due  to  the  absence  of  negative  terms  In  DCs). 
Therefore  the  zeros  of  D(3j  will  remain  In  the 
left-half  3-plane  and  low  frequency  unstaCle 
modes  cannot  arise  during  activation. 

Function  programmablllty ; 

The  objective  of  this  research  was  to 
develop  a  device  that  Is  capable  of  realizing 
the    following   transfer   functions:      LP  where  Ks) 

-  K/DCsi,    BP  where  T(3)    -  KS/D(3),    HP   where   T{s) 

-  KS^/Dls).  AP  where  T(s)-[32-3(u>p/Qp)«Up'' 1/Dls) 
and  N  where  T(3)-( 32*uj2)/d(  3) .  By  optimizing 
the  design  of  the  filter,  It  was  found  that  all 
of  the  above  functions  can  be  realized  by  the 
second  order  QIC  section  using  four  resistors, 
two  capacitors  and  two  OAs  as  shown  In  Table  1. 
These  passive  elements  are  connected  to  differ- 
ent nodes  to  achieve  the  various  realizations. 
A  set  of  CMOS  bilateral  switches  controlled  by 
digital  binary  word,  are  used  to  relocate  the 
same  elements  In  different  ways  to  achieve  the 
desired  filtering  functions  according  to  Fig.  2. 
The  truth  table  of  the  switches  control  logic  13 
shown  in  Table  2,  while  Fig.  3  Illustrates  the 
corresponding  minimized  CMOS  logic  circuit  used 
for   passive  elements  relocation. 

Paramater   programmablllty: 

While  four  of  the  resistors  are  equal  and  of 
value  R  each,  the  fifth  resistor  is  the  Qp 
determining  resistor  and  of  value  Rq-RQp.  The 
two  capacitors  are  equal  and  of  value  c-l/CupR) 
each.  Two  equal  banks  of  capacitors  are  used  to 
control'up.  Each  bank  contains  n  binary- 
weighted  capacitors  connected  In  series  with  n 
CMOS  switches  as  shown  In  Fig.  « .  Using  n  bit 
binary  word  to  control  the  switches,  2" 
different  values  of  o  can  be  obtained  that  cor- 
responds to  2"  different  valjes  of  up.  Using  a 
similar  technique,  the  value  of  Rq  can  be  con- 
trolled by  an  m  bit  digital  word  that  yields' 2™ 
different  values  of  Qp,  as  Illustrated  In  Fig.  5. 
Thus,  full  Independent  control  of  the  pole  pair 
ujp  and  Qp  are  achieved  by  programming  the 
digital  words  controlling  the  switches  to  obtain 
the  corresponding  c  and  Rq.  The  complete  second 
order  programmable  filter  Is  shown  In  F 1  g  .6.f. 
wnara  th*  function  prografflmabl 1 1 ty  as  well  •• 
the   parameter   programmablllty   are  demonstrated. 

Higher   order   programmablllty: 

Active  filters  design  procedure  can  be 
classified  as  direct  or  cascade.  In  direct 
synthesis  procedures  the  transfer  function  Is 
realized  as  a  single  section  [7].  In  cascade 
synthesis  procedures  a  high  order  transfer 
function  Is  expressed  as  a  product  of  first  and 
second  order  transfer  functions  and  each  of 
these  Is  realized  Independently.  The  overall 
network  Is  obtained  by  cascading  the  Individual 
sections.  The  cascade  method  of  synthesis 
offers  two  practical  advantages  (a)  simple 
network  tuning  (b)  a  few  number  of  universal 
sections  can  be  designed  which  can  realize  a 
multitude  of  network  specifications. 


The  second  order  CIC  network  structure 
lends  Itself  to  the  cascade  synthesis  procedure 
since  It  does  not  require  additional  Isolating 
amplifiers.  Flg.6.k  shows  a  block  diagram  of  a 
programmable  higher  order  filter  that  utilizes 
the  second  procedure  by  cascading  2  or  more 
sections  of  the  filter  network  shown  In  Fig.  6.3. 
The  result  Is  a  high  order  fully  programmable 
general  purpose  filter,  that  can  be  tailored  to 
match  almost  any   proposed  specification. 
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COMPUTER    SIMULATIONS     AND    EXPERIMENTAL 
VERIFICATIONS: 


Fig.  7  Shows  dlfferenot  computer  simulation 
outputs  of  the  programmable  filter.  The  plots 
simulate  the  filter  responses  assuming  ideal  OAs 
with  infinite  Cain  Bandwidth  Products  (GBWP),  as 
well  as  practical  filter  responses  assuming  OA's 
finite  GBHPs  of  1  MHz  as  of  that  of  the  LM7H1  OA. 
A  single  pole  OA  mode)  was  utilized  to  approxi- 
mate the  filter  transfer  functions  In  the  later 
case.  The  approximation  was  found  adequate 
since  the  simulation  results  of  the  nonldeal 
response  were  found  to  be  of  close  proximity  to 
the  experimental  results  of  Fig.  8.  The  experi- 
mental results  were  obtained  using  a  three  bit 
word  for  filter  topology  programmablllty  to 
select  the  type  of  transfer  functions.  A  two 
words,  four  bits  each,  were  used  for  filter  para- 
meters programmablllty  where  up  and  Qp  are 
controlled  Independently  as  given  In  Table  3. 
Fig.  7  also  Illustrates  a  higher  order  progran- 
mabllity  where  a  fourth  order  characteristics 
are  shown  for  a  LPF  and  a  Chebychev   BPF. 

1.      CONCLUSION 

The  novel  design  described  here  has 
resulted  In  a  universal  programmable  filter  that 
can  be  digitally  controlled  to  realize  almost 
any  practical  filter  specifications.  This  is 
done  through  the  use  of  CMOS  switches  controlled 
by  binary  codes  to  program  the  order  of  the 
filter,  the  filter  topology,  the  filter  center 
frequency  and  selectivity.  The  design  procedure 
required  developing  optimum  switching  arrange- 
ments for  the  minimum  redundancy  In  components 
and  the  least  dependence  of  the  filtering 
function  on  swltoMlng  Impar f co 1 1  ons  auoh  ■• 
switches  stray  capacitances  and  non-zero  and 
non-llnear  switch-on  resistance.  Further  invest- 
igation is  being  conducted  to  develop  a  program- 
mable switched  capacitor  realization  that  can 
allow  frequency  scaling  by  changing  clock 
frequency.  Work  Is  also  In  progress  for 
developing  an  extended  bandwidth  programmable 
filter  using  the  composite  operational  amplifier 
technique  proposed  earlier  by  the  author.  Such 
Implementation  would  lead  to  a  very  useful 
monolithic  device  at  moderate  cost. 
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Fig.  3  The  logic  controlling  figure  2  switches 
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